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Abstract — A LiNbCfe waveguide qua si-phase- matched 
sum-frequency generation (SFG) device using a fan-out domain- 
inverted grating for picosecond optical sampling was designed 
and fabricated. Sampling of signal waves at L5-^m wavelength 
by 25-ps sampling pulses generated by a gain-switched laser 
diode is demonstrated in a device of 5 mm length with a 
normalized SFG efficiency as high as 43%/W. 

Index Terms — Nonlinear optics, optical sampling, optical sig- 
nal processing, optical waveguide, quasi-phase matching, sum- 
frequency generation, wavelength conversion. 

I. Introduction 

NONLINEAR waveguide quasi-phase-matching (QPM) 
structures can perform wavelength conversion as well as 
various high-speed optical signal processing. Autocorrelation 
measurement of picosecond optica! pulses [1] and self phase 
modulation [2] in QPM waveguides have been demonstrated. 
Optical switching [3] and pulse compression [4] have been 
demonstrated using bulk QPM structures. 

High-speed optical sampling [5] using ultrashort optical 
pulses is an important technique for optical waveform mea- 
surements and demultiplexing of TDM optical signals. Ultra- 
fast sampling based on Type-Il phase matched sum-frequency 
generation (SFG) in bulk K.TP and organic crystals has been 
demonstrated [6], [7]. The SFG efficiencies, however, were in 
the order of 10" 2 %/W or less, and the experiments required 
optical amplifiers for signal and sampling pulses. Waveguide 
QPM-SFG devices would offer advantages of much higher 
efficiencies and larger flexibility for signal and sampling 
wavelengths. In this letter, we present design, fabrication and 
experimental demonstration of a LiNb03 waveguide QPM- 
SFG device for efficient picosecond optical sampling. 

IT. Device Description and Design 

Fig. 1 illustrates a LiNb03 waveguide QPM-SFG device 
for optical sampling. The device consists of an array of 
channel waveguide and a ferroelectric-domain inverted grating 
of a fan-out pattern [8]. A signal wave of frequency u?i 
(wavelength Ai) and sampling pulses of w 2 (7^1) (A2) are 
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Fig. 1. Schematic illustration of a waveguide QPM-SFG device for optical 
sampling. 
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Fig. 2. Dependence of the QPM grating period on the input wavelengths, 

combined together and coupled into a waveguide channel. The 
dependence of the QPM grating period on Ai, A2 in 1.5-fj.m 
band calculated by using the extraordinary index dispersion 
of bulk crystal is shown in Fig. 2 with A A = Ai — A2 as 
a parameter. At or near the QPM condition, sum-frequency 
(W3 = w\ fu/2) wave (A 3 ) of a power P3 proportional, to 
the product of the signal and sampling wave powers P u P 2 
is generated, and therefore the device serves as an optical 
multiplier for sampling. The SF wave can easily be detected 
by a near-infrared detector with appropriate time constant to 
give a sampling signal. 

For continuous-wave (CW) input waves, the F3 is given 
approximately by 

P 3 = n 2 P l P 2 L 2 {sin(AL)/AL} 2 J 

2A = 2n(N 3 /X z - N1/X1 - N2/X2 - 1/A) (1) 

where L is the device length, A the grating period, n the SFG 
coupling coefficient, Ni f 2 f z the effective indexes for Ai f 2,3 



1041-f I35/99$tO.OO © 1999 IEEE 



1028 



IEEE PHOTONICS TECHNOLOGY LETTERS, VOL. II, NO. S, AUGUST 1999 



qT 

f 



SE 




100 



4 8 8 

z Device length L [mmj 

Fig* 3. Calculated dependence of the normalized SFG efficiency and mini- 
mum pulsewidth on the device length. 



waves, and A denotes the deviation from exact QPM. The 
QPM bandwidth for SFG in terms of the grating period A is 
0.32, 0.064 for L = 1, 5 mm, respectively, as shown in 
Fig. 2, At the QPM-SFG condition (A = 0), weak (phase- 
mismatched) second-harmonic generation (SHG) of Ai, A 2 
waves may take place, and the SH waves affect the sampling 
performance as noise. The Ax/2 and A 2 /2 SH powers, Pu and 
Pzi* are given by formulae similar to (1) with SHG coupling 
coefficient k S kg (« = 2«shg for A! ~ A2) [9]. The extinction 
ratio may be defined as the ratio of residual SHG maximum 
(envelope) power to the SFG power (Piunv + Pii*uS)jP$* A 
large extinction can be obtained by choosing an appropriate 
wavelength separation AA depending on L. 

For SFG of short pulses, wavelength bandwidths of the 
pulses, pulse broadening in the dispersive waveguide, and 
walkoff between the input and SF pulses must be consid- 
ered Simple calculation assuming Fourier-transform limited 
picosecond pulses shows that, in a LiNbC>3 waveguide of 
L < 10 mm, the broadening and walkoff are so small that 
their effects can be neglected. For efficient SFG, however, the 
wavelength bandwidth must be within the QPM width which 
decreases with increasing L. Thus, the QPM width dominates 
the lower limit for the sampling pulsewidth. Fig, 3 shows the 
dependence of the normalized SFG efficiency at exact QPM, 
P3/P1P2 = K 2 L 2 t and the minimum pulsewidth upon A 
LiNb03 waveguide of 192 ^m 2 effective cross section 5^ [9] 
with a grating of 1:1 duty ratio is assumed. The result indicates 
that picosecond sampling with efficiency higher than 100%/W 
(~Q dB with l-W peak sampling pulses) can be obtained. For 
example, 63%/W SFG efficiency, 1-ps minimum pulsewidth 
and 1 8-dB extinction are expected for L = 5 mm, A A = 20 
nm and P1/P2 = 0.1. 

III. Fabrication 

A SFG device consisting of a fan-out grating of A = 
16.0-18.7 fim and 80 waveguide channels of L = 5 mm was 
fabricated. The fabrication process is similar to that of SHG 
and difference-frequency generation (DFG) devices reported 
previously [10], [11]. At first domain-inverted gratings were 
fabricated in a z-cut LiNb0 3 crystal of 0.15 mm thickness 
by applying a voltage pulse of ~4 kV between a corrugation 
electrode on the +z surface and a uniform electrode on the 
— z surface. The pulse duration was automatically controlled to 
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Fig. 5. Experimental setup for optical sampling. 

give a predetermined inversion charge. The grating has 1 : 1 
duty ratio and it continues through the entire crystal thickness. 
Then the waveguide array was fabricated by proton exchange 
in benzoic acid at 200 °C for 1.5 h using an Al film mask 
with 3.5-pm-wide channel openings followed by annealing at 
350 °C for 4 h in oxygen atmosphere. Finally, the waveguide 
input/output ends were polished for end-fire coupling. 

IV. Experimental Result 

The SFG characteristics were first measured by CW opera- 
tion. An external-cavity tunable LD and a DFB LD were used 
as light sources of Ai and A 2 , Both waves were combined 
by a 3-dB fiber coupler and coupled into the waveguide 
to excite fundamental TM modes. A channel was selected 
to accomplish QPM, and for A is A 2 = 1573, 1542 nm, 
strongest SFG of A3 — 779 nm was obtained in a channel 
of 16.47-^m grating period. The mode sizes (FWHM) were 
measured as 7.2 x 6.3 pm 2 for A x , A 2 and 2.4 x L8 fj,m 2 for 
A3 (£ e ff = 192 M m2 )- Fi g- 4 shows the dependence of P3 upon 
Pi, Pi measured at the waveguide output end. A normalized 
SFG efficiency of 43%AV was obtained. This value compares 
fairly well with the theoretically estimated value of 63%/W. 
The SFG bandwidth in terms of Ai was 4 nm, which was in 
good agreement with the theoretical value of 4.2 nm. 

Fig. 5 shows the setup for demonstration of optical sam- 
pling. A DFB LD was driven by a ^1 GHz sinusoidal 
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Fig. 6. Obtained sampling signal. 

wave with small dc bias to generate sampling pulses by gain 
switching. The pulses ware amplified by an Er-doped fiber 
amplifier. The average power of the coupled sampling pulses 
was roughly 10 mW. The pulsewidth was measured as 25 ps 
by a SHG autocorrelator. The wavelength and the bandwidth 
were A 2 = 1543 and 0.6 nrn, respectively. The external 
cavity LD was directly modulated by a 1-GHz sinusoidal 
wave to produce a signal wave of Ai = 1577 nm and 5-mW 
average power. The SFG output was detected by a Si p-i-n 
photodiode (insensitive in the 1.5-pm band) and displayed on 
an oscilloscope of 1-Mft input impedance. When the repetition 
frequencies of the signal wave and sampling pulses were set at 
1 GHz and 1 GHz-1 kHz, respectively, a sampling signal of 1- 
kHz repetition, which is a time-converted replica of the signal 
waveform, was obtained as shown in Fig. 6. The waveform 
directly observed by a high-speed InGaAs p-i-n photodiode 
and a sampling oscilloscope of input impedance is also 
shown for comparison. It should be noted that the level of the 
optical sampling signal was orders of magnitude higher than 
that of the direct detection. Observations of several different 
waveforms showed that the time resolution was consistent with 
the 25-ps sampling pulsewidth. The residual SHG level was 
smaller than 5% of the SFG sampling signal level. 
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V. Conclusion 

High-efficiency picosecond optical sampling using a 
LiNbCb waveguide QPM-SFG device has been demonstrated 
for the first time. Obtained efficiency was three orders of 
magnitude higher than that of sampling by bulk crystals. 
Experimental work is being continued for improvements of 
efficiency and time resolution and integration with couplers. 
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L5-//m-Band Wavelength Conversion 

Based on Cascaded Second-Order 
Nonlinearity in LiNb03 Waveguides 
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Abstract — We report wavelength conversion and spectral in- 
version using cascaded second-order nonlinearity in periodically 
poled L1N0O3 waveguides pumped at 1,5 /xm. The converter 
has an internal conversion efficiency of —8 dB, a conversion 
bandwidth of 76 nm, and a constant conversion efficiency for 
the 50-dB range of signal powers tested. 

Index Terms — Nonlinear optics, optical fiber communications, 
optical frequency conversion, wavelength conversion, wavelength- 
division multiplexing, ^ 

x 

WAVELENGTH conversion is an important function 
in wavelengt-division-|nultiplexed (WDM) optical net- 
works. Among numerous demonstrated wavelength conversion 
techniques, difference frequericy generation (DFG) [1] is at- 
tractive in several respects: it is an instantaneous process 
that easily accommodates mqre than terahertz modulation 
. handwidths, can simultaneously^ up and down convert multiple 
channels with equal efficiencies, has negligible spontaneous 
emission noise and has no intrinsic frequency chirp. 

DFG-based wavelength converters have been demonstrated 
in AIGaAs [1], [2] and LiNb0 3 waveguides [3]-[5], showing 
promising results for WDM wavelength conversion. However, 
the use of DFG and current device efficiencies indictate the 
need for a single mode pump laser with —50-100 mW of 
power operating in the 750-800-hm range. Moreover, it is 
difficult to simultaneously launch the 780-nm pump and 1.5- 
band signals into the fundamental mode of the waveguide, 
although this has been solved using an integrated mode 
coupling structure [4] In this letter, we demonstrate wavelength 
conversion in LiNb0 3 waveguides using a cascaded second- 
order nonlinearity (x t2) : x^ 2) ) [6]-[8], where both the input 
pump and signals are in the 1.5-fim band. Mode matching in 
this case is simplified compared to DFG pumped at 780 nm, 
since all input wavelengths are in the same band and can be 
launched into a single port. 

The ► x (2) -based device for 1.5-pm band wavelength 
conversion uses a pump in the 1550-nm region. The pump 
at frequency w p is upconverted to frequency 2w p by second- 
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harmonic generation (SHG) via the second-order nonlinearity 
X^. The generated 2w p simultaneously mixes with input 
signals ur t to generate wavelength shifted outputs a; ou t = 
2cj p — u a by DFG via another x^ process, Phasematching 
between interacting waves for both SHG and DFG is required, 
and this can be accomplished by choosing an appropriate 
quasi-phasematching (QPM) grating period. It is interesting to 
note that this x^ : X^ process mimics four-wave mixing 
(FWM) which uses the third-order nonlinearity x C3) - The 
effective x^ of such process in LtNb03 under quasiphase- 
matching condition is 10 4 -10 5 times larger than that of silica 
glass. The x^ : X^ allows the use of a very short sample 
when compared to fiber, and has better noise figure compared 
to FWM in semiconductor optical amplifiers. The converted 
output is related to pump power and input signal power in the 
low conversion efficiency approximation by [7] 

Pout * \v 2 L*P*P, (1) 

where Pv n %,Ppt and P a are converted output power, pump 
power and signal power, respectively, tj is the normalized 
efficiency (the same for SHG and^near degenerate DFG) in 
units of mW" l -cm" 3 (this number is proportional to both 
the mode overlap between the interacting waves and the 
material nonlinearity x^)> anc * *h e interaction length. 
The conversion efficiency scales with the fourth power of 
the interaction length, so doubling the device length will 
increase the conversion efficiency by a factor of 16. In the 
high-conversion efficiency regime, the exponent of L in 1 is 
less than 4 due to pump depletion. In an optimized LiNbC>3 
waveguide with a normalized efficiency of 150%/W-cm 2 (or 
0.0015 mW" l 'Cra" 3 ) and 6-cm interaction length, one can 
achieve 0-dB conversion efficiency with ^75 mW of pump 
power; or 3-dB gain with ~10G mW of pump power. 

We fabricated the waveguides by annealed proton exchange 
in periodically poled LiNb0 3 (PPLN) [9]. The device used 
in this experiment is 4 cm long, has a QPM period of 15 
Aim, waveguide width of 12 /im, proton exchange depth of 
0.7 fim, and was annealed for 26 h at 325 °C, The above 
parameters allow phasematching at room temperature between 
the fundamental mode of the pump at 1556 nm and the 
fundamental mode of SHG wave at 778 nm. We can tune the 
pump wavelength by using waveguides with different QPM 
period and/or temperature tuning. The normalized efficiency 
of this device is 65%/W-cm 2 . One can achieve a normalized 
efficiency of 150%/W-cm 2 by optimizing the mode overlap 
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Fig. 1. Experimental setup. ECL: External cavity laser. BPF: Bandpass filter. 
PPLN WG: Periodically poled LiNbOa waveguide. OSA: Optical spectrum 
analyzer. The waveguide is fiber coupled at both the input and the output. 
The later is analyzed by an OSA and power meter. 



between the fundamental mode of the pump and the first high- 
order mode of the second harmonic. At the input and output of 
this device,, the waveguides are optimized for fiber coupling by 
tapering a section of 1-rnm length to adiabatically transform 
the modes into and out of the wavelength conversion sections. 

Fig. 1 shows a schematic diagram of the experimental setup 
used in this work. The pump laser is an external cavity laser 
(ECL) amplified by an erbium-doped fiber amplifier (EDFA) 
to a level of —300 mW and filtered through a FBG in order to 
supress the amplified spontaneous emission (ASE). This pump 
is combined with signal sources generated in four different 
ECL/s and fiber launched into the waveguides. The output of 
the waveguide is fiber coupled and later analyzed in an optical 
spectrum analyzer (OSA) and power meter. The fiber to fiber 
coupling loss in this configuration is —4.5 dB due to reflection 
losses at the uncoated endfaces (—1.7 dB), mode mismatching 
between the fibers and the waveguide (—1.3 dB), and intrinsic 
waveguide losses (—1.5 dB, i.e., —0.35 dB/cm). 

We first characterized the device by measuring the SHG 
power versus pump wavelength. We chose to keep the wave- 
guide at 90 °C (or higher) in order to avoid photorefractive 
effects. This effectively shifts the phasematching wavelength 
to 1562 nm, but keeps the other parameters unchanged. The 
device shows a near-ideal sine 2 wavelength tuning curve with 
a peak internal efficiency (output SHG power divided by 
square of input pump power) of — 500%/W and an FWHM of 
—0.27 nm at low pump power. At higher pump power (>100 
mW), the device displayed a distorted wavelength tuning 
curve due to photorefractive effects of the generated second- 
harmonic waves, which results in waveguide nonuniformity 
and reduces the effective device length. 

We performed simultaneous multichannel wavelength con- 
version with a pump power of —110 mW launched inside the 
waveguide; the results, in Fig. 2, show that the efficiencies 
are the same (-15 dB) for all the input channels. The best 
conversion efficiency (-8 dB) was obtained with a pump 
power of —175 mW in a similar device with a slightly different 
QPM period and operated at 120 °C, as shown in the inset of 
Fig. 2. 

For : x^-based wavelength converters, the converted 
electric field is the complex conjugate of signal electric field, 
meaning that the output electric field spectrum is the mirror 
image of input spectrum about the pump wavelength. We show 
the spectral inversion properties of our device in Fig. 3 using 
an asymmetric input spectrum, formed by combining an ECL 
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Fig. 2. Simultaneous wavelength conversion of four WDM channels with 
the same conversion efficiencies. Pump power is ~1 10 mW inside the PPLN 
waveguide at a wavelength of 1562 nm. The inset shows a — 8-<fB conversion 
efficiency of a similar device operated at 120*0 with 175 mW of pump 
power. 
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Fig. 3. Spectral inversion of the '■ X* 2) wavelength converter. The 
output electric field spectrum is the mirror image of input spectrum about the 
pump wavelength, a feature that can be used to invert the signal chirp. 

signal and filtered ASE from /an EDFA. This feature can be 
used to invert the signal chirp for dispersion management in 
a transmission system. 

Fig. 4 shows the measured bandwidth of this wavelength 
converter when tuning the input signal wavelength while 
keeping the pump wavelength fixed at 1562 nm. This device 
has a 3-dB signal bandwidth of 76 nm, which is wider than 
the theoretical bandwidth of an ideal 4-cm-long device due to 
photorefractive effects at high pump powers. 

We tested the linearity of this wavelength converter by 
varying the input signal power, as shown in Fig. 5. We 
measure a linear response for more than 50 dB with a 
maximum input signal power of —0 dBm. The only saturation 
mechanism in this device comes from the depletion of SHG or 
pump power. We estimate a 0.16-dB deviation from linearity 
for a input power of 0 dBm in an ideal 6-cm-long device with 
a normalized efficiency of 150%/W-cm 2 and pump power of 
^100 mW. For a device with a 2.5-cm interaction length, a 
normalized efficiency of 65%/W-cm 2 , and pumped by —100 
mW, the deviation at 0-dBm signal power should be 0.006 dB, 
which is smaller than the resolution in this experiment. Our 
measurements confirm these estimates. 

The conversion efficiency of this device currently is limited 
by photorefractive effects at high pump powers, which can 
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Fig. 4. Conversion efficiency vs. input signal wavelength. The closed circles 
are the measured results and the solid line is a theoretical fit to a device with 
an effective interaction length of 2.5 cm. The 1562-nm pump has a power of 
~135 mW inside the waveguide. 




Input (dSm) 

Fig. 5, Single-channel output/input transfer curve of the x^ ' X^ wave- 
length converter. This device has a linear response for more than 50-dB range 
of input signal powers. 

be reduced by operating at higher temperature or annealing 
samples in an O2 atmosphere. In addition, it has been shown 
that MgO:LiNbt)3 has several orders of magnitude more 
resistance to photorefractive effects than congruent LiNbOs, 
and thus a similar device fabricated on MgO:LiNbC>3 would 
operate more efficiently. 

In summary, we have demonstrated wavelength conversion 
within the l.5-{im band using x C2) : X (2) in PPLN waveguides. 



This approach requires only lasers operating in the 1550-nm 
region. We have shown spectral inversion and simultaneous 
conversion of four WDM channels, a static internal conversion 
efficiency of -8 dB, a conversion bandwidth of 76 nm, and 
a linearity over an input power range of more than 50 dB. 
With further improvements, a wavelength converter with 0- 
dB conversion efficiency or even gain is a feasible goal. In 
addition, this '• -based wavelength converter can also 
function as clock recovery or time gating component by use 
of a pulsed pump. 
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All-Optical Switching Based on Cascading of 
Second-Order Nonlinearities in a Periodically Poled 
Titanium-Diffused Lithium Niobate Waveguide 

Hirohisa Kanbara, Hiroki Itoh, Member, IEEE, Masaki Asobe, Kazuto Noguchi, Member, IEEE, 
Hiroshi Miyazawa, Tsutomu Yanagawa, and Itaru Yokohama 



Abstract — We report efficient all-optical switching by using a 
periodically poled titanium-diffused lithium niobate (Ti:PPLN) 
waveguide. The periodically domain inversion of the Ti:PPLN 
waveguide is achieved by electric field poling. The switching 
is driven by the second-order nonlinear effect of the cascad- 
ing of phase-matched sum-frequency-generation and difference- 
frequency-generation processes. It is found that about 12% of the 
signal is switched by a 20- W gate power. 

Index Terms — All-optical switching, cascading, difference- 
frequency-generation, periodically poled lithium niobate, 
second-order nonlinearity, sum-frequency-generation. 



I. INTRODUCTION 

PERIODICALLY POLED lithium niobate (PPLN), which 
provides large second-order nonlinearity and high trans- 
parency, has led to the development of optical devices that are 
based on the high-yields frequency mixing. Typical examples 
are blue light sources using second-harmonic-generation and 
communication-band- wave length converters using difference- 
frequency-generation (DFG) [1], [2], Besides these devices, 
PPLN is also promising for constructing efficient optically 
gated optical switches [3], It has been pointed out that the 
cascading of the second-order nonlinearities makes possible 
a lower gate power compared with that for the conventional 
switch driven by third-order nonlinear effects like a nonlinear 
change in the refractive index [4]-[6]. As an instance of 
this, we proposed optical switching that utilizes the cascading 
of phase-matched sum-frequency-generation (SFG) and DFG, 
and thereby realized full-switching operation with a gate power 
of 3 kW using 10-mm-long bulk PPLN [3]. In this study, for 
further improvement of the switching efficiency, we fabricated 
a titanium-diffused PPLN (Ti:PPLN) waveguide and examined 
the optical switching properties based on the cascading of 
phase-matched SFG and DFG in the waveguide. 
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Fig. 1. Preparation process for the Ti:PPLN waveguide. The periodically 
domain inversion results from electric field poling. We used a 300-^m-thick 
substrate for the poling. 



II . Experimental 

The Ti:PPLN waveguide was prepared by using electric field 
poling [7], [8], The preparation process is shown in Fig. 1. A 
titanium-diffused lithium niobate (LN) waveguide was formed 
on the -z face of the LN substrate [9]. For the poling, an 
aluminum electrode grating was deposited on the +z face of 
the substrate, and the insulator was also coated with photoresist 
on this +z face so that its domain would remain unchanged. 
The period of the electrode grating was 15.5 fjxti. A lithium 
chloride liquid electrode directly touched the — z face of the 
substrate [7]. After the application of electric field, we found 
that the domain inversion only occurred in the part where the 
aluminum grating was deposited. A 10-mm-long waveguide 
was used to conduct the optical switching experiment. 

In optical switching hy means of the cascading of SFG and 
DFG processes, SFG occurs first, when a strong gate beam, 
^Gate, and a weak signal beam, ^signal, are input together, 
and then DFG yields the switched signal 



t^Gate + ^Signa.1 — * W 5FG 
WSFG — k>Gate — * ^Signal* 



(1) 

(2) 
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Fig. 3. Signal wavelength dependence of the SFO efficiency. The signal 
with a wavelength of 1.556 provides the largest SFG efficiency for a 
1.319-/*m gate wavelength. 

The experimental setup is illustrated in Fig- 2. We examined a 
three-terminal switching operation. The configuration is almost 
the same as that of an optical Kerr shutter switch, which 
employs a third-order nonlinear process. When there is no 
gate beam, the signal does not pass through the analyzer. 
On the other hand, when the gate beam strikes the Ti:PPLN 
waveguide, the signal passes through the analyzer because the 
direction of the polarization of the emerging phase-changed 
signal is effectively rotated [3]. In case the phase change of 
the signal is 7r, a 1 00% switching efficiency is obtained. We 
constructed a collinear-type switch with a dichroic mirror. For 
the gate beam, we used a fiber-compressed, mode-locked YAG 
laser pulse with a wavelength of 1.319 pm. The gate pulses 
had a duration of 6 ps and a repetition rate of 100 MHz. A 
CW distributed-feedback laser diode coupled to an erbium- 
doped fiber amplifier (EDFA) produced the signal beam. The 
largest second-order nonlinear component, d^, of LN was 
parallel to the z axis. For the switching experiment, the gate 
beam was linearly polarized parallel to this z axis, and the 
polarization of the linearly polarized signal beam is set to 
7r/4 with respect to that of the gate beam. The optical filter 
prevented the gate beam from entering the germanium p-i-n 
photodetector, which was connected to a lock-in amplifier. The 
waveplate compensated for the birefringence of the waveguide. 

III. Results and Discussion 

Before trying any switching, we first did the SFG mea- 
surement to determine the appropriate signal wavelength. In 
the SFG measurement, the signal beam was produced by a 
wavelength-tunable semiconductor CW light source and an 
EDFA. In Fig. 3, the SFG efficiency is plotted as a function 
of the signal wavelength. The dots are measured data, and 
they show that the SFG efficiency strongly depends on the 




"0 TO 20 30 
Gate power (W) 

Fig. 4. The dependence of the switching efficiency on the gate peak power. 
The signal was found to continuously increase with increasing gate power. 




t/4 0 Jt/4 n/2 3™*4 
Gate poliirization angle, 6 

Fig. 5. The dependence of the switching efficiency on the angle of polariza- 
tion of the gate beam. The z-directional gate beam gives the largest switching 
efficiency, but no signal was detected for the y-directional gate beam. 

signal wavelength. The main peak was found to be at 1.556 
ftm, so we used that wavelength in the switching experiments. 
The solid curve was calculated from the Sellmeier coefficient 
at 300 K [10], but the center wavelength was shifted so 
that it matched the experimental data. The refractive index 
difference between LN and Ti:PPLN probably brings about 
this discrepancy in the center wavelength. Large error bars 
for experimental data are thought to be due to the instability 
of the output power of the wavelength-tunable semiconductor 
CW light source. The measurement revealed that the maxi- 
mum normalized SFG conversion efficiency for the Ti:PPLN 
waveguide was 130%/W. This value is more than two orders 
of magnitude larger than the 0.97%/W for bulk PPLN [3]. 

The switching property is shown in Fig. 4. This graph shows 
how the intensity of a switched signal varies with the gate peak 
power. The gate power was estimated at the output of the 
waveguide. The lens-to-lens transmittance of the waveguide 
was 1 5% for the gate beam and 20% for the signal beam. The 
extinction ratio of the linear polarization maintenance for the 
signal beam was over 30 dB. The signal power was 1 mW at 
the output of the waveguide. We observed clear switching of 
the signal even when the gate power was less than 20 W. At the 
gate power of 20 W, the switched signal was approximated to 
be 12°/o. This result indicates that the gate power is over two 
orders of magnitudes smaller than that for the conventional 
optical Kerr shutter switch using a silica glass fiber [11]. To 
obtain the same switching efficiency, the needed gate power 
was found to be 5.5 times larger than the calculated one. This 
large power consumption is probabEy caused by nonuniformity 
of the periodically poled structure, optical damage of the 
waveguide, and incomplete overlapping of the gate and signal 
beams. Trials for improvement are in progress. 

Fig. 5 shows how the switching efficiency depended on the 
angle of polarization of the linearly polarized gate beam. The 



330 



inset shows the waveguide and the direction of polarization of 
the gate beam. Here, B = 0 means that the gate polarization is 
parallel to the z axis, or in other words, parallel to the^as of 
LN, The polarization of the signal beam was kept at B - tt/4. 
We obtained the maximum switched signal when the gate beam 
was polarized parallel to the z axis. And we did not get any 
signal when the gate beam was polarized in the y direction. For 
the waveguide we used, the wavelength conversion efficiency 
is much lower in the y direction, because there is a phase 
mismatch in this direction and the second-order nonhneanty 
of d-ji is small Assuming zero nonlinearity along the y axis, 
the calculated curve agrees well with the measured data. This 
proves that the switching was driven by cascading processes 
of the second-order nonlinearity governed by d 3 3 as well as 
phase matching, and that the switching did not originate from 
the third-order nonlinearity where phase matching is irrelevant. 

IV. CONCLUSION 
All-optical switching that are based on the cascading of 
phase-matched SFG and DFG was demonstrated by use of a 
10-mm-long TuPPLN waveguide. We constructed an optical- 
Kerr-shutter-type configuration to perform a switching experi- 
ment and observed an efficient three-terminal optical switching 
operation. The switched signal was about 12% at the gate 
power of 20 W, It was evidenced that a PPLN waveguide will 
be very useful in producing highly efficient optical switching 
devices. 
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Low-Power All-Optical Gate Based on Sum 
Frequency Mixing in APE Waveguides in PPLN 

KL R. Parameswaran, M. Fujimura, M. H. Chou, and M M Fejer 



Abstract — We present an all-optical gate implemented in pe- 
riodically poled lithium niobate. Efficient mixing is achieved by 
using a phase-matched guided-wave interaction. A control wave 
at 1.537 Mm is used to gate a signal at 1.552 fim, where a control 
power of 185 mW is sufficient to achieve 96% depletion of a low- 
power signal. A simple switch configuration is described whereby 
high-contrast low-power all-optical switching can be performed. 

Index Terms — Gated mixer, nonlinear optics, optical fiber com- 
munications, optical frequency conversion, optical switching, op- 
tical waveguide, periodically poled lithium niobate, quasi-phase- 
matching, sum frequency generation. 



I. Introduction 

ALL-OPTICAL switching is an enabling function for fu- 
ture high-speed fiber communication systems. Previously 
demonstrated approaches using the third-order material non- 
linearity x (3) suffer from difficulties such as high-switching 
powers and/or long devices (due to weak nonlinearities in op- 
tical fiber-based four wave mixing [1]) or signal degradation 
due to additive noise (in four wave mixing in semiconductor 
optical amplifiers (SOA's) [2]. Cascading of second-order non- 
linearities (x (2) : X^ 2) ) has also been used t0 perform switching 
[3]. The switching powerneeded in these configurations is quite 
high (on the order of several Watts), primarily because they op- 
erate far from phase matching in order to simultaneously ob- 
tain a large phase shift (required for switching) and flat spec- 
tral response. 

We demonstrate an all-optical gate based on sum frequency 
mixing (SFM) using a single x C?) interaction in an annealed 
proton exchanged (APE) waveguide formed in periodically 
poled lithium niobate (PPLN). This approach exploits the 
high conversion efficiency available from a phase-matched 
interaction, resulting in 96% depletion of a CW signal by a 
185-mW control beam. A structure using this gate is described 
with which high-contrast low-power, all-optical switching can 
be performed. 
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EL THEORY 

Coherent three-wave mixing using x (2) has several attractive 
features such as transparency to signal format and the addition 
of negligible excess noise. Difference frequency mixing (DFM) 
using has been used to perform efficient wavelength 
conversion within the 1.55-^m communication band [4] as 
well as between the 1.3- and 1.55-pm bands [5], In each case, a 
strong pump wave is mixed with weak signal waves to produce 
mixed outputs at frequencies mirrored about half the local 
oscillator frequency. These same devices can be used to mix 
a low-power signal with a stronger control wave to generate 
a signal at the sum frequency. Both phenomena are described 
by the well-known coupled mode equations for three-wave 
mixing [6], but with different boundary conditions in each 
case. In DFM, the low-power signal is amplified during the 
mixing process, whereas in SFM, the signal wave is depleted 
during generation of the sum frequency wave. In the simple 
case of a phase-matched interaction without loss, where the 
control wave is much stronger than the signal wave (and can 
be considered undepleted), the evolution of power in the SFM 
and signal waves is described by 

Fsfm W = Aig(O) sin 2 (v'w^CTRLl') (la) 



Fsig(L) = ftio(O) cos 2 (Vw-Pctrl£) Ob) 

where L is the interaction length, Fctrl is the control power, 
A sfm and A S cg are the SFM and signal wavelengths and 77 nor is 
the normalized conversion efficiency. Complete depletion of the 
signal occurs when the argument of the trigonometric functions 
equals tt/2. Hence, gating of the signal can be accomplished by 
turning the control power on and off. 

III. Experiment 

The device used in this work was designed to mix signal and 
control waves in the 1.55-^m band, with a mixing region de- 
signed to be insensitive to variations in waveguide width, which 
loosens fabrication tolerances [7]. This noncritical design re- 
quires a waveguide that supports multiple transverse modes at 
both input wavelengths in the conversion region, so a mode filter 
and adiabatic taper at the input are incorporated to facilitate 
launching of power into the fundamental modes in the conver- 
sion region [8J. The optimized widths are 4 p,m for the mode 
filter and 12 £im for the conversion region. The 6.05-cm long 
sample was diced from a 3-in-diameter wafer of LiNb03, which 
had been electric-field poled [9] with a quasi-phasematching 
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Fig. 2. Typical pulse traces, showing depletion of the signal in the presence of 
the control pulse. 



(QPM) period of 14.7 pm. The QPM grating is 5.55 cm long. 
Waveguides were formed by proton exchange for 15 h at a tem- 
perature of 160 °C (to a depth of 0.71 fxm), followed by an- 
nealing for 26 h at 328 °C. The end faces of the sample were 
then polished to permit efficient end fire coupling. 

The experimental setup is shown in Fig. 1. Control pulses 
were produced by a fiber-amplified externally modulated diode 
laser. A DFB laser diode operating at 1.552 ^m produces 
the CW signal. Phasematching was observed with Actrl 
= 1.537 ^m resulting in Asfm - 0.772 ^m. The signal and 
control beams are combined, then launched into the same 
waveguide. Output light at the three wavelengths is separated 
using a dichroic mirror and diffraction grating, then directed 
onto fast photodiodes for power measurement. The experiment 
was carried out at room temperature. 

IV. Results and DrscussioN 

Fig. 2 shows depletion of the CW signal in the presence of the 
control pulse. One advantage of the nonlinear transfer function 
of the SFM process is that distortions in the control pulse are 
suppressed in the depleted signal pulse, making the switching 
less sensitive to variations in the control pulse. Fig. 3 shows 
a plot of normalized signal transmission as a function of con- 
trol power. The solid line shows the calculated result obtained 




50 100 150- 200 
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250 



Fig. 3. Measured and calculated variation of signal transmission with control 
power. 

by numerical integration of the coupled mode equations de- 
scribing the SFM process in waveguides, including the prop- 
agation losses neglected in obtaining (1). The calculation used 
a normalized internal SFM conversion efficiency of 33.9 W -1 
(as determined by measurement of the second harmonic gener- 
ation efficiency of 8.47 W^ 1 ) and typical values of propagation 
losses at the three wavelengths (0.35 dB/cm at \$ig and Actrl > 
and 0.70 dB/cm at Asfm) - Mode profiles used in the calculation 
were obtained by solution of Maxwell's equations subject to the 
refractive index profile resulting from the annealed proton ex- 
change process [16]. Nearly complete (96%) extinction of the 
signal is seen at an input control power of 185 mW, a value in 
reasonable agreement with theory. The mechanism responsible 
for the residual 4% is unclear; slight phase mismatch arising 
from waveguide nonuniformity is one possibility. 

V. Future Possibilities 

The required gating power in SFM devices can be signifi- 
cantly reduced by several relatively straightforward means. The 
device used in this experiment had nonuniformities resulting 
from imperfect fabrication. Optimization of the process would 
result in a 30% increase in 7? nor . The conversion efficiency in 
waveguide QFM interactions is roughly proportional to the 
square of the grating length. Periodically poled devices as long 
as 8 cm have been reported [10], hence, moving to this length 
from the current value of 5.55 cm would increase 77 nor by about 
a factor of two. High index cladding layers [U] and buried 
waveguides [17] have also been used to improve the spatial 
overlap of modes in waveguides. Calculations indicate that this 
technique could potentially increase 7} nol by another factor of 
two. Combining the above techniques could result in a control 
power around 40 mW (or 16 dBm, a power level readily avail- 
able from commercial EDFA's). This corresponds to a gating 
energy of 1 pJ for 25 ps pulses in a 20-GHz return-to-zero pulse 
train, which compares very well with similar experiments using 
SOA's [12]. The bandwidth of the device is limited by group 
velocity walkoff between the three interacting waves such 
that the minimum pulsewidth scales inversely with the device 
length. Hence, efficiency and bandwidth can be traded against 
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Fig. 4. "AND" gate using x ca> in NOLM configuration. 
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each other. One way to eliminate this tradeoff is known as 
"quasi-group velocity matching" [8], where after each walkofT 
length, the faster pulse passes through a delay line integrated 
into the waveguide structure in order to be resynchronized with 
the slower pulse. This scheme, currently being investigated, 
would allow for high efficiency switching of "arbitrarily" short 
pulses. 

VI. Proposed Switch Structure Using x (2) Gate 

Higher contrast can be achieved by placing this device in 
an interferometer or in a nonlinear optical loop mirror (NOLM 
[13], Fig. 4). In the absence of the control pulse, the loop is 
balanced and all of the signal power exits at port 1. When the 
control is on, signal power traveling in one direction in the loop 
gets depleted so that incomplete interference at the 3-dB cou- 
pler results in signal power leaving port 2 (at a level 6 dB below 
the input value). Hence, an "AND" logic function is implemented 
(there is output only when the signal and control are both on). 
Power switched to the SFM wave (at 0.772 pm in the current 
device) can be used for further signal processing, A similar de- 
vice has been used to perform efficient optical sampling [14]. 

An attractive feature of SFM is that the signal wavelength ac- 
ceptance bandwidth is narrow (—0.3 nm for a 6 cm long device) 
such that individual WDM channels can be switched out of a bit 
stream in this way. In order to switch multiple wavelengths si- 
multaneously, DFM can be used as the mixing process instead of 
SFM. Here, the NOLM is unbalanced by amplifying the signal 
rather than switching it out Parametric gain on the order of 3 dB 
has been observed in APE waveguides in PPLN [15]. This can 
be used to create a switch with 0-dB insertion loss and very wide 
signal bandwidth (around 46 nm, for a 6-cm device). 

VII. Conclusion 

An optical gate based on phasematched SFM in QPM-APE 
waveguides in PPLN has been demonstrated. Nearly complete 
extinction with a control power of 185 mW was observed. Two 
possible configurations for very high contrast switching have 
been presented. 
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proton exchanged Lithium Niobate waveguides. In z-cut crystals, the 
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1. Introduction 

Among man-made crystals for electro-optics and nonlinear photonics, Lithium Niobate (LN) is 
probably the most investigated and employed. It allows efficient phase and amplitude 
modulation and parametric generation in guided-wave configurations [I]. Several specific 
technologies have been mastered in order to fabricate LN planar integrated devices, from 
Titanium in-diffusion to ion-exchange to e-beam implantation, etc. Among them, proton 
exchange (PE) combines the advantages of a low-temperature process with a substantial 
surface index increase and a reduced photorefractive sensitivity [2]. Despite their high quality, 
even with additional improvements such as thermal annealing, soft exchange in diluted 
solutions, sealed ampoule method, etc., the asymmetric profile of PE waveguides is an 
obstacle to efficient coupling to/from silica fibers. In order to improve such coupling between 
waveguide and LPfjl fiber modes, Reverse PE (RPE) was recently introduced to obtain buried 
and more symmetric channel waveguides.[3-8] 

RPE is a two-step process: a first, deeply penetrating rr*->Li + exchange followed by a 
second, weaker Li + — >H* exchange. Such technique allows to form two superimposed 
waveguides with graded index profiles, each of them able to confine one of the principal 
(orthogonal) polarizations. For a z-cut crystal, for instance, an ordinary ray can be confined as 
a TE guided mode, while an extraordinary ray becomes a TM eigenfields. As sketched in 
Fig. 1, TE and TM waves can then propagate and carry power along y at different depths, 
corresponding to the two (ordinary and extraordinary) refractive distributions. 



2.20 2.25 2.30. 
Of— r " 




i \ TE-guide(FE) ,L 



TM-guide (SH) 



substrate 



Fig. I. RPE geometry with graphs of ordinary and extraordinary index distributions in z-cut LR 

These superimposed waveguides can be coupled by a nonlinear interaction involving a 
non-diagonal element, I e., thru a polarization field orthogonally polarized with respect to the 
input electric wave. This is the typical configuration of a Type I phase-matched process of 
Second-Harmonic Generation (SH-G), whereby the harmonic is orthogonal to the fundamental 
frequency (FF) field. Furthermore, this is the standard SHG taking place in uniform 
birerringent phase-matched LN crystals thru the d3j element of the quadratic susceptibility, the 
largest non-diagonal. It is worth underlining that the double exchange tends to recover the 
quadratic nonlinear response of the surface layer, contrary to single-step PE [9]- Therefore, 
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using mode coupling between fundamental and second harmonics in the surface and buried 
guides of an RPE geometry, respectively, one expects the generation of a TM frequency- 
doubled field by injecting a TE FF wave in the upper waveguide. 

In this Paper we investigate numerically SHG in RPE planar waveguides in z-cut 
propagation-homogeneous (non-periodically-poled) LN, and demonstrate that large conversion 
efficiencies can be obtained even in those cases when much higher or lower temperatures 
would be required in the absence of Quasi-Phase-Matching, with good fabrication tolerances 
and physical separation of the harmonic field distributions at the output. 

2. Models 

Let us consider a z-cut y-propagating LN crystal. The pertinent non-diagonal elements of the 
second-order tensor d will couple z and x components of the electric-field. An FF wave 
propagating in the surface waveguide as a TE mode will radiate a cross-polarized SH field into 
the buried guide, supporting TM eigenwaves (see the extraordinary -dashed- index profile in 
Fig. 1). We will assume the upper guide is single-moded at the FF wavelength, whereas the 
deep guide is multimodal at the SH, the most common situation in a practical realization: 

E? E (z,y) = A(y)f,(z)e^ N ^ , E&(z, y) = Z v B v (y)f 2 v .(z)e- i2k ° N - y (1) 

with ko the vacuum wave-number at FF and N the effective index. The generation will then 
take place through the interaction between the fundamental FF mode and high-order SH 
modes, whenever phase-matching is accomplished with an appreciable overlap integral 
I = \ (f 2 v J*f *dz , as prescribed by coupled-mode theory (CMT) [10]. 

In a structure of finite length, however, the interaction is better described by a system of 
partial differential equations taking into account the details of the geometry and the overall 
field distributions undergoing nonlinear propagation and energy exchange. For the 
superposition of cross-polarized (FF and SH) harmonics propagating along y in a z-cut LN 
crystal, using Maxwell equations in the usual approximation of weak guidance and keeping the 
second derivatives with respect to y, we get: 

+ 3^ + k^(z)) 2 E^ E + 2M J5 = 0 

£5™ + ^£*™ + 4k^"(z)) 2 E^ + 4kk, (E?e> 2 - 0 

dy (r4r> dz 

with and n 2 ™ the ordinary and extraordinary graded index distributions at FF and SH, 
respectively, and njj* the bulk (i. e., substrate) ordinary and extraordinary indices at SH, 
d 15 =d 3l the relevant nonlinear coefficients when Kieinmann symmetry holds. The set (2) can 
be integrated using the Enhanced Finite Difference Beam Propagation Method (EFDBPM) 
with transmitting boundary conditions [1 1]. 

To evaluate the nonlinear response in a realistic case, we considered typical index profiles 
for planar RPE waveguides, as prepared and reported by Korkishko et al. in z-cut LiNb0 3 [8], 
The samples were first PE-treated for 14h in an Ammonium Di-Hydrophospate melt at 220°C, 
then they were thermally annealed for 1 1 Oh at 330°C; finally » they were RPE-processed in a 
L1NO3 melt for lOOh at 300°C. For a wavelength X=27tc/co=1.32jim such that the top guide 
supports one mode at FF, i. e., the index distributions can be expressed as: 
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n£(z)= n^+AnSexpfja.)* (3 ) 

with q=l (2) for the FF (SH), i=l (2) V 2<2 0 (Vz>^), z 0 =6^im, G 0 , =4um f a o2 -5u.rr;, 

a el =2um, a c2 =3|xm, A% = - 0,023 and An^= 0.1. 

For the profiles (3) and q=l, we calculated the FF eigenmodes with a Finite Difference 
Method (FDM). Then, we integrated eqns. (2) launching in the upper guide suitable FF field 
distributions and studying the y-evolutions of FF and SH waves. In our simulations we 
employed either an LPOl-like (x-in variant) Gaussian beam with a symmetric distribution 

u^z) = exp[-(z-z,) 2 fo 2 } optimally centered at z^uxn with oH).83z,; or the zeroth- 
order TE eigenmode as derived numerically by FDM, Integrating the SH field and its intensity 
over the output transverse (z) section at a distance L from the input, we evaluated SHG 
conversion efficiencies (per unit wavefront along x) for various excitation levels, temperatures 
and device lengths. 

3. Results and discussion 

Let us first compare the outcomes from standard CMT and EFDBPM with eqns, (2). For a 
I cm-long sample at 25°C, it is apparent from Fig, 2 that the computed conversion efficiency 
with the two approaches (using a Gaussian excitation for the EFDBPM) is markedly different. 
Examining in more detail the case of a Gaussian excitation, mapping FF and SH intensities 




Input power (W/'um) Input power (W7u.m) 

Fig. 2. Computed SHG conversion efficiency versus input power for a 1cm sample at 25°C (a) 
CMT and (b) EFDBPM results In case of Gaussian excitation. 

versus propagation and depth, we see in Fig. 3 that, for a 10W/pn input power into an RPE 
waveguide of length L=2cm, nearly 80% of the input beam couples to the TEq mode of the 
upper waveguide, with the remaining 20% exciting a leaky or quasi-mode, as apparent from 
the lobes visible across the depth coordinate z (Fig. 3a). The parametric nonlinearity generates 
a TM polarized SH field which, confined in the buried waveguide, grows with propagation. 
The competition between the two excited modes at FF, however, limits the overall up- 
conversion efficiency. This insight is supported by Fig. 4, showing the effective index versus 
modal number at SH at 25°C. The 10 mode is close to phase match the fundamental FF 
mode, but the FF quasi-mode can interact via higher-order modes, thereby limiting the amount 
of frequency doubled power through Maker-like interference. Moreover, due to the particular 
geometry and in contrast to the collinear guided-wave case, the overlap integral with the TEo 
at CO is an increasing function of the modal order at 2co, as displayed in the inset of Fig, 4, In 
any case, the EFDBPM clearly provides a more accurate estimate of the up-converted energy 
at the output. 
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y(cm) y( c "X> 

Fig. 3, Contour maps of (a) FF and (b) SH intensity versus z and y for a Gaussian input at 
lOW/jim in a. sample at temperature 25 < *C 




Modal order 

Fig, 4. Phase-matching diagram at 25°C versus modal order m at SH. The horizontal tines refer 
to the FF modes, TEo and leaky (or quasi-mode), respectively. The inset shows the 
corresponding TM ra -TEo overlap integral. 

The conversion process can be enhanced by exciting only one FF mode, launching an input 
beam with transverse profile mapping the TE Q mode. In this case, at an optimum temperature 
of 85 °C corresponding to phase-matching, the SH grows monotonically along y largely 
depleting the fundamental, as visible in Fig. 5. 

Using this optimum input profile, by a transverse integration in y=L we computed SHG 
conversion versus input FF power (in W/urn). The results are shown in Fig. 6 for two different 
sample lengths and temperatures. Due to the nature of the interaction and the substantial 
pump depletion, the conversion efficiency is sub-linear with both FF excitation and square of 
the propagation distance. For L=2cm it reaches 90% at lOW/um, with a maximum normalized 
value of 14% |im/Wcm (Fig. 6b, solid line). The temperature dependence, taken into account 
at each wavelength and polarization [12], is rather weak. This is clearly visible in Fig. 7, 
where the SHG efficiency is graphed versus temperature for two different input powers and 
sample lengths. 




y(cm) y( cm ) 
Fig. 5. Contour maps of (a) FF and (b) SH intensity versus z and y for a TEo. input at lOW/fim. 
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Input power (W/pjn) Input power (W /jjjti) 

Fig. 6. Conversion efficiency versus input FF power (TE<, excitation) at (a) 25 and (b) 85°C, for 
samples 1cm (dashed line) and 2cm (solid line) in length. 




87 J J 50 as r»,3 150 



Temperature (°CT) Temperature (°C) 

Fig. 7. Conversion efficiency versus temperature (TEg excitation) at input powers of (a) 
1 W/um and (b) 10W/|im, for samples 1cm (dashed line) and 2cm (solid line) in length. 



4. Conclusions 

The major drawbacks of Cerenkov and temperature-tuned SHG, i. e. conical emission and 
critical thickness in PE waveguides [13-14] or critical temperature/modal tuning in Thin- 
diffused channels [15], respectively, can be eliminated in propagation-uniform Reverse Proton 
Exchanged guides in z-cut Lithium Niobate. The conversion efficiencies, numerically 
estimated using a propagator based on a robust PDE model, are as high as 90% at a 
fundamental wavelength of 1.32um in 2cm-long planar guides optimally excited at lQW/\im. 
Such values compare well to the best results reported to date in non-periodically-poled LN (L 
e,, non quasi-phase-matched) waveguides, both in Cerenkov [16-17] and guided- wave 
interactions [15, 18]. 

RPE waveguides appear to be promising candidates for robust and low temperature 
frequency doubling, with the unique advantages of SH confinement, weak temperature 
sensitivity, polarization selection and physical separation of the harmonics at the output. We 
are currently investigating channel waveguides for bidimensional confinement and 
longitudinally tapered transverse profiles for improved coupling/handling of the up-converted 
output distribution. 
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We report depth profiling of the <f 33 nonlinear coefficient in annealed proton exchanged LiNbC>3 
waveguides using reflected second-harmonic generation from angle-lapped samples. At depths 
greater than the initial proton exchange depth the </ 33 coefficient retains nearly its bulk LiNbOj 
value, but within the initial proton exchange region the value of d 2 $ varies with annealing time. 
These results reconcile previous conflicting measurements of the d$ 3 coefficient and explain the 
variation in the efficiency of guided wave frequency conversion devices with annealing. 



Quasi-phase matched (QPM) frequency conversion in 
LiNb0 3 waveguides is an efficient method for the genera- 
tion of short wavelength radiation from infrared laser di- 
odes. While different techniques have been used to form 
the ferroelectric domain grating necessary for QPM, the 
annealed proton exchange (APE) process has been use4 
exclusively to provide waveguide confinement. Contradic- 
tory measurements of the nonlinear optical properties of 
proton exchange (PE) and APE-LiNb0 3 waveguides have 
been previously reported. 1 " 6 The d 33 coefficient of PE- 
LiNb0 3 has been measured to be between 0 (Refs. 4 and 
5) and 0.63 (Ref. 3) of the bulk LiNb0 3 value. The effect 
of annealing is similarly controversial, with measurements 
of d 33 in APE-LiNb0 3 Varying from 0% (Ref. 4) to 90% 
(Ref. 3) of the bulk value. Recently demonstrated 7 * 8 
guided-wave QPM second-harmonic generation (SHG) 
devices with normalized conversion efficiencies > 600%/ 
W cm 2 indicate that APE-LiNb0 3 waveguides must have a 
large nonlinear* coefficient. 

In this letter we report measurements of the depth 
dependence of the d^ nonlinear coefficient in APE-LiNb0 3 
waveguides using reflected SHG from angle-lapped sam- 
ples. With a 532 nm fundamental wavelength, the 266 hm 
second-harmonic (SH) wavelength is above the APE- 
LiNb0 3 band edge so that only the SH generated within a 
0*05 iitsx absorption depth 4 of the surface is observed. The 
samples were angle lapped so that lateral position could be 
mapped onto depth into the sample, allowing depth profil- 
ing via lateral translation. This technique circumvented 
complications caused by the graded refractive index profile 
and the unknown spatial variation in the J 33 coefficient in 
the APE-LiNb0 3 . The depth profiles presented here ex- 
plain the variation in the d iZ values reported in the litera- 
ture. Normalized conversion efficiencies calculated with 
our results are consistent with those observed from guided- 
wave QPM SHG devices. 

APE waveguides were fabricated on x-cut LLNb0 3 by 
proton exchange in pure benzoic acid at 173 D C for 66 min 
to a depth of 0.42 fim, calculated using the results of Ref. 
9, and annealed in air at 333 °C for varying times. The 
samples were then polished at a wedge angle < 2 mrad, 
with the actual relationship between lateral position and 
depth determined using surface profilometry. 532 nm radi- 
ation from an injection seeded 2-switched frequency- 



doubled Nd-YAG laser was focused onto the sample at 
normal incidence with a 30 /im FWHM spot and a peak 
intensity of <150 MW/cm 2 . The. samples were scanned 
under the focused spot using a motorized translation stage 
and the reflected 266 nm SH was detected using a dichroic 
mirror, a solid blind photomultiplier, and a gated integra- 
tor* The fundamental and SH fields were polarized parallel 
to each other and the z axis of LiNb0 3 and thus were 
coupled by the <tf 33 nonlinear coefficient. To discriminate 
against polishing artifacts, each sample had an internal, 
unexchanged reference portion formed by masking half the 
sample with Al prior to proton exchange. The mask was 
removed by etching in NaOH prior to annealing. A portion 
of the sample remained unwedged, so that reflected SH 
could be obtained from the original surface. SHG and pro- 
filometry measurements were performed over the whole 
surface to eliminate errors associated with polishing skew. 
Figure 1 shows the sample orientation and geometry, with 
the dashed line indicating the beginning of the wedged 
portion of the sample and the origin for measurements of 
lateral position; negative values of lateral position repre- _ 
sent data from the unwedged portion of the sample, i.e. t 
the original surface. 

Figure 2 shows the reflected SH power at 266 nm, 
normalized to the signal from bulk LiNb0 3 versus depth 
and lateral position for APEWLiNb0 3 waveguides annealed 
for 3, 5, 9, and 63 h. The error bars represent the variation 
in the depth measurement for 5 different scans. The figure 
inset shows the spatial step response of the detection sys- 
tem projected onto depth, measured by scanning the fo- 
cused spot off the end of a bulk LiNb0 3 sample. This 




FIG. 1, Geometry of the wedged sample used for reflected 
SHG measurements. 
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FIG. 2. Reflected 266 nm SH intensity, normalized to that from bulk 
IiNbOj/vs lateral position (top scale) or depth into the waveguide (bot- 
tom scale) for samples proton exchanged to 0.42 pm and annealed for 3, 
5, 9, and 63 h. The inset shows the step response of the 
detection system. 

demonstrates a depth resolution of <0.05 /Ami The samples 
displayed constant reflected SH over the entire unex- 
changed side with intensities equal to that from bulk 
IiNb0 3 /with the exception of the sample annealed for 63 
h, which had a polishing artifact between lateral positions 
—75 and +50 /zm. 

For the APE-LiNb0 3 sample annealed for 3 h no re- 
flected SH signal was observed from the original surface or 
from the portion of the wedged surface corresponding to 
the original PE layer. However, there is an abrupt increase 
from 09& to 80% of the bulk LiNb0 3 value at a position 
corresponding to the interface between the original PE film 
and the LiNb0 3 substrate. At lateral positions correspond- 
ing to depths greater than the initial PE depth, the reflected 
SH power approaches that measured from bulk LiNbC 3 . 
At depths >1 fim the measured reflected SH power is in- 
distinguishable from that of bulk LiNb0 3 . 

With further annealing, the position of the abrupt in- 
crease in the reflected SH power moves toward the original 
surface, and the reflected power approaches that from bulk 
LiNfa0 3 at shallower positions. No reflected signal was 
ever observed from the original surface regardless of an- 
nealing time, except from the sample annealed" for 63 h at 
lateral positions where polishing damage was observed. 
However, no signal was observed over the rest of the un- 
wedged portion of this sample. 

The two previous measurements 4,6 of the cf 33 coefficient 
of APE-LiNb0 3 waveguides performed using reflected 
SHG with a 532 nm fundamental source were made on 
unwedged samples, only probing a 0.05-/£m-thick region at 
the original surface. Reference 2 reported no reflected SH 
regardless of annealing time at annealing temperatures of 
350 *C, similar to our results when the probe is confined to 
the unwedged part of the sample. Reference 4 reported that 
APE-£iNb0 3 proton exchanged for 0.5 h and annealed for 
1.5 h at 310 *C or for 1 h at 350 *C exhibited reflected SH 
with an intensity of 1/2 or 1/4 that obtained from bulk 



LiNb0 3 . For samples proton exchanged longer than 0,5 h, 
no reflected SH was observed, regardless of annealing time 
or temperature. 

The absence of any reflected SH from the unannealed 
PE region could be explained by the formation of a PE- 
LiNb0 3 film with either inversion symmetry or a different 
X m tensor that did not result in any reflected SH. Simi- 
larly, the progression of the abrupt stepiike recovery to the 
surface region with increasing annealing could be ex- 
plained by the solid state epitaxial regrowth of the IiNb0 3 
phase from the original interface between the PE film and 
the LiNb0 3 substrate. After the abrupt increase in the re- 
flected SH from 0% to 80% of the bulk value there is a 
slower approach to the bulk value occurring over the next 
0.5 ftm in depth. In this region the measured d 33 value 
correlated with concentration, independent of annealing 
time, so we conclude that for low concentrations the *f 33 
coefficient is slightly reduced from its bulk value. While the 
details of the progression of the g? 33 recovery to the original 
surface may be dependent on the orientation of the LiNb0 3 
substrate and annealing temperature, it is unlikely that the 
*f 33 coeflicient of APE-LiNb0 3 waveguides fabricated on z- 
and jc-cut substrates differ substantially at depths greater 
than the initial PE depth. Understanding the absence of 
reflected SH from the original surface of APE-LiNb0 3 
may require microstructural and x-ray diffraction experi- 
ments to elucidate the crystal structure within the top 0.05 

of the surface. This is of little practical importance 
since the modal amplitudes arc very small at the surface. 
To design guided wave frequency conversion devices it is 
adequate to assume that at depths greater than the initial 
PE depth the cf 33 .coefficient nearly retains its bulk value 
throughout^ annealing, while variations occur within the 
original PE depth. 

The measurements of the d 33 coeflicient shown in Fig. 
2 can be used to determine the normalized conversion ef- 
ficiencies (77) for any guided wave nonlinear optical inter- 
action- As aa example, we compute 77 for SHG of ^=850 
nm radiation in 2-cut, xr-propagating APE-LiNb0 3 
waveguides fabricated like those shown jn Fig. 2. We as- 
sume first-order QPM over an interaction length / with a 
uniform, depth-independent ferroelectric domain grating, 
which can be fabricated using an electric-field 7 or electron- 
beam* poling technique. With the output SH power given 
by Pto^vtfjP* V k given by 

where rf QPM = 2d^°Vir, dg*(z) = d$iHz)/d\i*°\ 
and the fields are normalized to carry unity power. The 
depth dependence of the refractive index profile was deter- 
mined using a one-dimensional nonlinear diffusion model 
for the APE process, 10 and we assume a 4-^im-wide top-hat 
dependence for the lateral refractive index profile. We as- 
sume that the two-dimensional modes are separable, 
E(yj)=E y (y)E z (z), and compute the mode profiles at 
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FIG. 3. Fundamental a„=850 nm) and SH (A^=425 am) modes 
(dotted lines), modal overlap (solid line), mealed normalized nonlin- 
ear coefficient 2? 3 re , the modal overlap with rfJf E (shaded region), and 
the H + concentration profile vs depth In APE-LiNhOj wayegwdes an- 
nealed for (a) 3 h and (b) 9 h. The abscissa is the same for (a) and (to, 

A =850 nm and A^=425 nm using the effective index 
method. The depth" overlap integral was evaluated using 
the depth dependence of the </ 33 coefficient derived by tak- 
ing the square root of the normalized reflected SH power 
given in Fig. 2. Shown in Fig. 3 versus depth are the optical 
modes, the measured normalized nonlinear coefficient 
2$? E (z), the modal overlap £(s)^»C*>. ^ as the 
shaded region, the modal overlap with the measured value 
of the nonlinear coefficient, d$% E (z)E i <a iz)E 2o (z), for 
waveguides annealed for 3 and 9 h; the data for a 5 h 
anneal are not shown for. brevity. Also shown is the H 
concentration profile, proportional to the index profile. At 
& the surface index changes are 0.0503, 0.0364, and 
0.0230, for 3, 5, and 9 h anneals, respectively. The overlap 
using the measured value of d^ yields n versus annealing 
of 410, 950, and 1090%/W cm 2 . For comparison, the 
modal overlap assuming the d 2 coefficient to be indepen- 
dent of depth yields r} versus annealing of 2010, 1640, and 
1280%/Wcm 2 . 11 For unannealed waveguides 17 =0 Since 



rf,,=0 for unannealed PE-LiNb0 3 . For anneal times 
longer than 9 h the modal overlap is removed from the 
region where <f 33 =0 and 17 decreases because of reduced 
confinement. Figure 3 shows that the variation m 17 with 
annealing is not due to a restoration in the <f 33 coefficient of 
the APE-LiNbOj waveguide as much as a variation in the 
overlap integral due to a spatial redistribution of the modal 
fields away from the initial PE region. . 

There are two recent reports of guided-wave QPM 
SHG using depth-independent ferroelectric domain grat- 
ings, with 7j of 600 (Ref.7) and 700%/W cm 2 (Ref. 8) for 
devices with APE waveguides fabricated similarly to the 
ones modeled above. These values fall within the range 
modeled, indicating that the general form of the d 33 depth 
profile presented here allows reasonably accurate device 
modeling. However, waveguide loss, laser longitudinal 
mode structure, QPM grating and waveguide inhomogene- 
ities and the value of the. J 33 coefficient of bulk LiNbOj . 
(recently measured to be 23.7 pm/V 12 rather than the ac- 
cepted value of 34.5 pm/V) render the modeling of the 
absolute i\ of real devices difficult and questionable. Mea- 
surements of 17 as a function of waveguide processing are 
ongoing to refine the models for the linear and nonlinear 
optical properties of APE-LiNb0 3 . 
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Abstract 

We report the first realisation of periodic domain inversion by electric-field poling in Ti:LiNb0 3 waveguides* for 
quasi-phase-matching (QPM) applications* Using a tunable Ti:AI 2 0 3 laser we demonstrate guided blue light generation in a 
third-order QPM interaction. 

Keywords: Lithium niobate: Quasi-phase-matching; Blue light generation 



1, Introduction 

Second harmonic generation (SHG) by quasi-phase- 
matching (QPM) in LiNb0 3 , via periodic modulation of 
the nonlinear coefficient, provides an attractive route for 
the realisation of blue and green light sources with several 
mW of output power* Electric-field poling has recently 
been used extensively for periodic domain reversal in 
ferroelectric materials, and hence modulation of the non- 
linear coefficient. This technique allows the use of the 
large d 53 coefficient and provides deep, laminar domain 
structures with high potential conversion efficiencies. Sev- 
eral devices have been demonstrated using this technique 
for QPM-SHG from IR wavelengths to the visible, both in 
bulk geometries [1,2], and in waveguide structures [3,4]. 
However, as with most other waveguide QPM-SHG de- 
vices demonstrated thus far [51, the waveguides in Refs. 
[3,4] were fabricated using the annealed proton exchange 
technique. The only demonstration of periodic domain 
inversion for QPM-SHG in a Ti:LiNb0 3 waveguide to 



date has been that by Webj5m et aL [6J, using a periodi- 
cally segmented waveguide. Domain inversion in this case 
was achieved through outdiffusion between the patterned 
Ti segments. However, the conversion efficiency obtained 
was 0.1 %/W- cm 2 , much smaller than that predicted by 
theory, possibly due to the domains being too shallow to 
allow a good overlap with the optical modes. Extension of 
electric-field poling to H-tndiffused waveguides would 
circumvent this problem by providing deep domains, hence 
allowing better interaction with the guided modes, and 
would also provide additional design options over those 
available with proton-exchanged waveguides, such -as al- 
lowing the use of both the extraordinary and ordinary 
waves for frequency conversion, and offering opportunities 
for active guides doped with rare-earths. In this paper we 
demonstrate, for the first time to our knowledge, blue light 
generation by QPM in an electric-field poled TiiLiNbO* 
channel waveguide. 
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2. Fabrication 

Channel waveguides were fabricated parallel to the 
X-axis by diffusing Ti stripes, 95 nm thick and ranging in 
width from 3 u,rn to 16 u.m in 1 u.m steps, into the — Z 
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Fig. I. Etched — Z face of a section of the sample, showing a 10 
|xm wide waveguide along the X-axis. The non-uniform domains 
within the waveguide are clearly visible. 



face of a 0.2 mm thick piece of LiNb0 3 (Z-cut) at 1005*0 
over a period of 9 hours. Our initial attempts to pole 
high-temperature-metal-diffused LiNb0 3 had been unsuc- 
cessful owing to the presence of a thin Li 2 0 outdiffusion- 
induced domain-inverted layer on the +Z-face t as revealed 
by wet etching, which blocked the e-fleld poling. In this 
case, to facilitate poling we removed the outdiffused layer 
by polishing off 50 |im of material from the H-Z face after 
the waveguide fabrication, thus ensuring that the crystal 
domain structure was unidirectional. Next, 1 mm . wide 
rectangular slots of period 9 j^m and a duty cycle of 33% 
(3 \Lva lincwidth), were opened in a photoresist mask on 
the — Z face of the sample, oriented perpendicular to the 
waveguides. This period was chosen as suitable for third- 
order QPM frequency-doubling to the blue from an IR 
wavelength of approximately 850 nm. Poling was achieved 
by applying four pulses of ~ 3 kV and *■* 300 ms duration 
via liquid electrodes, as described previously [lj. Fig. 1 
shows the -Z face of a small part of the sample, etched in 
a mixture of HF and HN0 3 . The photograph shows a 10 
u.m wide waveguide, with evidence of non-uniformity of 
the grating, including random fluctuations in. the domain 
wall positions, within the waveguide. Furthermore, the 
final grating structure obtained in the waveguides had, on 
average over the complete interaction length, a duty cycle 
between 60% and 80%. The occurence of this non-optimal 
mark/space ratio requires further investigation to allow 
accurate control of device parameters. A similar grating 
pattern was observed on the +Z face, the etched regions 
on this face corresponding to the non-etched regions on the 
—Z face, hence demonstrating continuous poling through 
the sample, with the rectangular domain shape generally 
maintained. The waveguide endfaces were then polished, 
yielding 5 mm long guides, periodically poled along the 
entire length. 



3. Optical characterisation 

A tunable Ti:Al 2 0 3 laser was used for the characterisa- 
tion of SHG from the waveguides. In this initial study we 
investigated blue light generation from a 6 u,m wide guide. 
Fig. 2 shows the third order SH output, which emerged 
entirely from the waveguide, as a function of the funda- 
mental wavelength. The fundamental light beam in this 
case was polarised along the z-direcn'on, Le. TM polarised, 
producing an interaction with the d 33 coefficient. As shown 
in the figure, conversion was achieved from the TMqq 
mode at the fundamental frequency to the TMqq and TM 10 
modes at the SH frequency. It was noticed that, for the 
same input power, the conversion efficiency to the second 
order harmonic mode was higher than that to the first 
order. We believe this to be due to a nonuniform! ty in the 
grating across the waveguide over a substantial part of the 
device length, as revealed through close inspection of the 
device, so that the second order mode interacts with the 
fundamental over a QPM structure which has a duty cycle 
closer to that required for optimum conversion efficiency 
in third order QPM. Moreover, there may be a reduced 
nonlinearity near the surface, due to the diffusion process, 
and variations in the index profile along the length of the 
waveguide, which would have a strong Detrimental effect 
on the coupling to the TM^ mode. Index variation effects 
are also evident in the difference in shape of 'the phase 
matching curves for the first and second order mode 
interactions. The phase matching curve for the TMJq — ► 
TMqq interaction is centred at — 845 nm, in good agree- 
ment with the calculated wavelength using dispersion rela- 
tions for TLLiNb0 3 waveguides [71 The curve however 
shows a FWHM bandwidth of —0.25 nm, twice that 




835.0 837.5 840.0 842.5 845.0 847.5 850.0 
Wave le ngth/n m 

Fig. 2. Guided second harmonic generated power normalised to 
the launched pump power, as a function of the fundamental 
wavelength. The phase matching wavelength in the bulk was 
found to be ~ $34 nra. 
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calculated from theory. This could be caused either by the 
irregular grating period or the above mentioned variation 
in effective index along the guide, due to inhomogeneities 
in the waveguide [8,9]. Fig. 3 shows the measured blue 
power in the fundamental mode of the harmonic frequency 
as a function of the coupled fundamental power at — 845 
nrn, the line in the figure being the best quadratic fit to the 
experimental points. The maximum power obtained in the 
blue was 7.9 uAV for a coupled power of 45 mW, giving a 
conversion efficiency of 1,6%/W « cm 2 . The conversion 
efficiency in this non-optimised demonstration device is 
more than one order of magnitude lower than that expected 
from theory. This, as mentioned above, is believed to be 
due to the grating mark/space ratio not being ideal, 
non-perfect axial and transverse uniformity of the grating, 
non-optimised overlap between the interacting modes and 
inhomogeneities in the effective index. 

The experimental points in Fig, 3 deviate slightly from 
the expected quadratic dependence, indicating the presence 
of photorefractive damage. Although periodic poling has 
been observed to increase the photorefractive damage 
threshold of LiNb0 3 [lJ. Goldberg et aL have recently 
observed damage effects at room temperature in their SHG 
experiments [2], These were attributed largely to a non-ideal 
grating structure, which would also apply to our wave- 
guide. At coupled fundamental powers greater than 45 
mW, the photorefractive damage was even more evident, 
with the blue output radiation spreading along the Z-axis. 
No attempt was made, in this initial experiment, to heat the 
sample and thereby anneal the damage. However, it is 
expected that it may be possible to reduce the photorefrac- 
tive damage in periodically poled H:LiNb0 3 waveguides 
by moving to a period smaller than the mode size in the 
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Fundamental launched power/mW 

Fig. 3, Second harmonic generated power as a function of the 
launched pump power. Note that the SHG power is that detected 
in the fundamental SHG mode at ~ 422.5 nm. The solid line is 
the best quadratic fit to the measured data points. 
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Fig. 4. Phase matching curve for the interaction of the fundamen- 
tal mode in the TE polarisation to the second order mode of the 
second harmonic in the TM polarisation (TEJo -» TM ). 



waveguide and ensuring that the domain structures are 
continuous and well defined, Using a period smaller than 
the mode size would enable a better compensation of the 
accumulated charges at the top and bottom of the beam 
profile [10], These charges have opposite signs in adjacent 
domains as the photogalvanic current, which produces this 
charge transport across the beam, follows the direction of 
the Z-axis. The resulting transverse electric field compo- 
nent, and therefore the induced refractive index change via 
the electro-optic effect, would thus be reduced. 

Finally, the polarisation of the input mode at the funda- 
mental wavelength was rotated by 90°, thus exciting the 
guides in the TE polarisation, and an interaction was made 
with the <tf 3l coefficient The phase matching curve for the 
TEqJ> -» TMjo interaction, centred at 876 nm, is shown in 
Fig. 4. The theoretically predicted wavelength for conver- 
sion to the first order SH mode is ~ 884 nm, which was 
outside the tuning range of our laser. The experiment 
however demonstrates the versatility of poled Ti:LiNb0 3 
waveguides in terms of accessibility to different wave- 
lengths using the same waveguide. Moreover, using the 
ordinary polarisation for a given order of QPM-SHG for a 
given blue wavelength would require a longer period of 
domain reversal than needed for the extraordinary polarisa- 
tion, alleviating* some of the problems encountered in 
fabricating short periods. Frequency" conversion using the 
d 31 coefficient is less efficient (d 3l ~ <U7<f 33 ) but may be 
important for interactions which involve very short pulses 
[II], as the group velocity mismatch between interacting 
fields at different frequencies is reduced in this polarisa- 
tion scheme, due to the lower dispersion. This makes the 
effective pulse interaction length longer, thus maintaining 
a good d cfS l e{f product. Possibilities for RE-doped wave- 
guide lasers in LiNb0 3 can also be extended, as the lasing 
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signal in these waveguides may be In the ordinary (<x) 
polarisation. 

4. Conclusions 

We have demonstrated, for the first time to our knowl- 
edge, periodic poling in a TuLiNbOj channel waveguide 
using electric-field techniques, and generated blue light via 
third order QPM. Further work is in progress to optimise 
the domain uniformity and modal overlaps within the 
waveguide, and achieve control over the resulting 
mark/space ratio. With these optimisations, and in particu- 
lar by using a first order QPM interaction, significant 
increases in conversion efficiencies are expected. More- 
over, we have recently demonstrated periodic poling in a 
Nd-diffused LiNb0 3 substrate [12] and also in a Tm-dif- 
fused Ti:LiNbO a waveguide [13], and intend to combine 
RE^diffusion and poling in Ti-indiffused waveguides for 
intra-cavity frequency doubling within a waveguide laser. 
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Fabrication, Characterization and Index Profile 
Modeling of High-Damage Resistance Zn-Diffused 
Waveguides in Congruent and MgO: Lithium Niobate 

W. M. Young, M. M Fejer, M. J. F, Digonnet, A. F. Marshall, and R. S. Feigelson 



Abstract — A study or the fabrication and optical properties 
of planar waveguides fabricated in MgO:LiNbOj and LiNb0 3 
substrates by diffusion of a ZnO film is presented. Transmission 
electron microscopy was used to show that using ZnO instead of 
metallic zinc as a source, and maintaining the ZnO film thickness 
below a prescribed value, greatly reduces second phase precip- 
itation and produces usable waveguides. Dopant and refractive 
index profiles were characterized by electron microprobe analysis 
and interference microscopy, respectively. The dependence of the 
Zn diffusion coefficient on temperature and the dependence of 
the refractive-index change on Zn concentration are inferred 
from these measurements. A simple model is also reported which 
predicts the Judex profile of the waveguide given the film thick- 
ness, diffusion time and temperature. The validity of the model 
is demonstrated by comparison between calculated profiles and 
profiles measured by prism coupling and IWKB analysis. 

I, Introduction 

OPTICAL waveguides in Lithium niobate (LiNb0 3 ) have 
been widely studied for applications including telecom- 
munication systems, nonlinear optics and fiber sensors. To date 
the two most widely used waveguide fabrication processes 
are titanium diffusion and proton-exchange (PE) [1], [2]. 
Both processes yield waveguides with low propagation loss, 
typically lower than 0.5 dB/cm [2], [3]. However, neither 
method produces waveguides with ideal performance for all 
applications. li-diffused waveguides in both LiNb0 3 and 
MgO:LiNb0 3 guide both polarizations but suffer from pho- 
torefractive damage even at wavelengths as long as 0.85 
^m [4], [5]. On the other hand, PE waveguides are much 
more resistant to photorefractive damage but guide only the 
extraordinary polarization. 

The objective of this work was to develop a new type of 
waveguide in LiNbOs which can guide both polarizations and 
exhibit high resistance to photorefractive damage. The current 
model for the photorefractive effect involves a refractive-index 
perturbation due to an interaction between the electrooptic 
effect and the local electric field created by charge separation 
under visible laser illumination [6]. Since the charge genera- 
tion and trapping centers are strongly related to dopant valence 
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state we speculated that cation dopants with a single valence 
state less than or equal to 2, such as Mg + * and H\ will have 
similar effects on photorefractive damage. This hypothesis was 
based on the experimental observations that both Mg* and 
H + reduce the photorefractive damage [7], [8], while cation 
dopants with multivalent states or a valence state greater than 
or equal to 3, such as Fe + W 2 , Nd +3 and IT 4 , increase the 
photorefractive damage [9], [10], 

In this work, Zn was chosen as a dopant because it matches 
the requirements predicted from the above arguments. In a 
recent article, we reported the fabrication and optical properties 
of planar waveguides made by the diffusion of ZnO into both 
LiNbOj and MgO:LiNb0 3 - Waveguides made in the latter 
material guide both polarizations, although with a substantial 
difference in index profile, and have propagation losses in the 
range of 0.3 to 1.2 dB/cm. They also exhibit no photorefractive 
in-plane scattering up to 90 kW/cm 2 at 514.5 nm, which is one 
of the highest values reported in the literature [11]. 

In this paper we present detailed material and optical studies 
of Zn-diffused waveguides in both MgO:LiNb0 3 and LiNb0 3 . 
Fabrication conditions required to obtain low loss waveguides 
are identified, in particular to avoid the formation of second 
phase precipitation of LiZnNb0 4 and surface degradation. We 
report measurements by electron microprobe analysis of the 
dopant distribution in Zn-diffused lithium niobate waveguides, 
as well as optical characterization of their index profiles. Thjs 
data is used to obtain values of the diffusion coefficient of Zn 
in congruent and MgO-doped LiNb0 3> their dependence on 
temperature, and the relationship between the refractive-index 
change and the Zn concentration. We also introduce a simple 
model, employing the diffusion and fabrication parameters, to 
calculate the dopant concentration profile, and from it the index 
profile of the waveguide. Comparison between the predictions 
of this model and the measured waveguide index profiles 
establishes the validity of this model. 

IL FABRICATION 

A. Experimental Procedures 

To prepare the Zn-diffused samples, a ZnO film was first 
sputtered onto the LiNb0 3 substrates. The dimensions of the 
substrates used to produce waveguides were 20 x 6 x 2 mm 
for MgO:LiNb0 3 and 20 x 6 x 1 mm for LiNb0 3 . Diffusion 
studies for both materials were carried out on 6 x 4 x 1 
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ram substrates. During diffusion, the samples were placed on 
platinum foil stands in the center of a covered alumina boat 
filled with congruent lithium niobate powder. The purpose of 
the powder is to suppress the net transport of Li ions in or out 
of the lithium niobate crystal [12]. To control the Li 2 0 vapor 
pressure in the vicinity of the LiNb0 3 , we used a mass of 
powder substantially larger than that of the sample (about 13 
g versus less than 0.6 g). The alumina boat was placed in the 
center of an alumina process tube, closed at both ends with 
alumina caps, which, in turn, was placed in a tube furnace 
for diffusion at elevated temperatures in dry air [13]. The 
temperature difference between the control thermocouple and 
a monitor thermocouple at the sample location was less than 
aS^C. In the text the quoted temperatures are provided by 
the control thermocouple. The thermal cycle, carried out in 
dry air, typically involved heating at about 8°Omin up to the 
diffusion temperature of 1000-1 100° C, diffusion for typically 
less than 1 hour, and cooling at an initial rate of lO°C/min. 



B. Surface Roughness 

To produce low scatter loss waveguides it is crucial that 
the surface quality be not degraded by the diffusion process. 
In our initial diffusion experiments the surfaces of all the 
samples were translucent after diffusion. Simitar observations 
of surface degradation resulting from the diffusion of metallic 
Zn into UNb0 3 were previously reported by Yoon et al [14 J. 
Based on the Zn0-Li20-Nb2O s pseudo-ternary phase diagram 
[15], we speculated that the surface roughening originated 
from the formation of a new phase on the substrate surface. 
The phase diagram suggesis that LiZnNbO* is a stable com- 
pound, which could therefore form near the sample surface and 
degrade its quality. For comparison, in the case of the diffusion 
of Ti in lithium niobate, the Ti0 2 -Li 2 0-Nb20 3 pseudo-ternary 
phase diagram indicates that T\Oi and LiNb0 3 can coexist in 
equilibrium al temperatures up to 1100°C [16]. Consequently, 
during the diffusion of Ti some metastable phases may form at 
elevated temperature but they will dissolve into lithium niobate 
over the course of the fabrication process [17], [18], 

Transmission electron microscopy (IBM) confirmed the 
presence of LiZnNbO* precipitates near the diffusion surface. 
Fig, 1 shows a bright field TEM picture of the substrate 
surface after diffusion of a 400 nm-thick film of ZnO at 
1000°C for 1 h. Isolated second phases, on the order of a 
few microns in size, were observed across the surface of the 
substrate. Energy-dispersive spectrometry confirmed that these 
precipitates contained high Zn concentration. The electron 
diffraction pattern of the precipitate has been identified as 
LiZnNbO*. Fig. 2 shows the electron diffraction pattern of 
the substrate (Fig, 2(a)) and of the precipitate (Fig. 2(b)). The 
lattice constants measured in Fig. 2(b) are consistent with doc- 
umented values for LiZnNbO*. The precipitates are always in 
a (201) orientation and have a specific in-plane alignment with 
respect to the (001) LiNbOj substrate. The {100} and {112} 
planes of the precipitates are parallel to the {10.0} planes of 
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Fig 1. Bright-field TEM photograph of a Zn-diffused congruent LiNbOj 
surface showing a precipitate (400-nm ZnO diffused at 1000° C for 1 hour). 



the substrate. For the LiZnNbO* structure, the angle between 
the two sets of planes is 59.2° , whereas the in-plane angles for 
<001> LiNbOj are 60°, The reflections from the precipitate 
(Fig. 2(b)) can be superimposed on the substrate (Fig. 2(a)) 
which indicates that the precipitates are slightly strained to 
give better atignment with the substrate. The orientation re- 
lationship between precipitates and substrate also gives three 
variants for the orientation of the precipitates and all of these 
can be observed within a given precipitate region. These results 
strongly suggest that a LiZnNb0 4 precipitate formed epitaxi- 
ally on the top surface of the waveguide during Zn diffusion. 
These isolated precipitates were probably responsible for the 
prohibitively high loss observed in early waveguides. 

There are at least two approaches to eliminate these second 
phase precipitates: 1) synthesize the dopant source from a com- 
position within the LiNb0 3 -2nNb 2 06-ZnLiNb0 4 tie triangle 
where LiNbOs is in equilibrium with quaternary Zn containing 
compounds or 2) find a set of parameters for ZnO diffusion 
which avoids the precipitation. The first approach would 
involve an in-depth study of dopant synthesis techniques and 
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Fig. 2. Electron diffraction patterns from the LiNbOj subsrraie (a) and 
prccipitarc (b) measured from the sample shown in Fig. }. The reflections 
indexed {30.0} planes from the subsirate LiNbOj in Fig. 2(a) matches the 
{400} and {224} planes from the precipitate in Fig. 2b, indicating epitaxial 
relationship between the substrate and the precipitate. 

compositional control during deposition of the source on the 
substrate. As a result the second alternative appeared to be 
a more desirable approach for establishing the usefulness of 
Zn-diffused LiNbOj waveguides. 

The key point in the second approach is to avoid any 
kinetic paths leading to compound formation. To this end 
ZnO was chosen instead of metallic zinc to avoid unnecessary 
metal oxidation which may enhance the formation of the 
intermediate phase [19], We also found that an important 
parameter controlling the formation of the precipitate is ZnO 
film thickness. We established experimentally that compound 
formation is greatly reduced, if not completely eliminated* 
when the film thickness is kept below some critical value. This 
critical thickness is approximately 160 nm for MgO:LiNb0 3 
and 100 nm for LiNb0 3 . No precipitate was observed in 
these materials with thicknesses of ZnO smaller than these 
values and for diffusion temperatures anywhere between 1000 
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and 1100°C. Satisfying the thickness condition was found 
to be essential to maintaining a good surface quality and 
producing low loss waveguides. By appropriate selection of 
the diffusion time and temperature (between 1000 and 1100°C 
for MgOrLiNbO* and around 1000*C for congruent UNbOj), 
the Zn concentration in the diffused region is high enough 
(typically 2-3 mol9£>) to provide good optical confinement, 
and yet low enough to maintain a good surface quality. 

III. Dopant And Diffusion Characterization 
Electron microprobe analysis (EMPA) to measure the dopant 
profiles in the Zn-diffused region required longer diffusion 
times (l-10h) and thicker ZnO films (0.4-1.0 m) than those 
used for the fabrication of waveguides. This was necessary 
so that high enough Zn concentrations would be present over 
greater depths to permit accurate EMPA measurement. Prior 
to characterization, a Zn-diffused sample was first sandwiched 
between two glass slides, then cut into two pieces. Each piece 
was polished at both ends, perpendicular to the diffusion 
surface. One piece was used for EMPA measurements, the 
other for Mach-Zehnder interference microscope measure- 
ments, as described in Section IV. The glass slides were 
used to prevent rounding of the edges of the diffused region 
during polishing, which would otherwise upset EMPA and 
microscopy measurements. 

An example of dopant profile is shown in Fig. 3 for an x- 
cut MgO:LiNbOa substrate after diffusion of a 200-nm ZnO 
film at 1000° C for 10 h. For all diffused samples studied 
here the ZnO film was not completely depleted at the end 
of the diffusion process. This assumption was checked by 
integrating the measured dopant profile along the depth of 
the diffused region to obtain the total amount of Zn present 
in the sample, and comparing it to the amount of Zn initially 
contained in the film. It is well known from diffusion theory 
that when the film is undepleted (and the diffusion coefficient 
is independent of concentration) the profile is a complementary 
error function (erfc). This is in agreement with measured 
profiles, as illustrated by theoretical fit of the experimental 
profile to an erfc plot in Fig. 3. Consequently, the diffusion 
coefficient D was obtained by fitting each measured dopant 
profile with an erfc function. 

Table I lists the values of D obtained from a total of 
10 samples in 3 different types of substrates (z- and x- 
cut MgO;LiNb0 3 and z-cut congruent LiNb0 3 ) processed at 
temperatures between 900 and HOFC. For short diffusion 
times, typically 1 hour or less, a nonnegligible amount of 
diffusion takes place during both the heating and cooling 
parts of the cycle. This was taken into account in correcting 
the measured diffusion coefficients by first estimating D(T), 
without applying any corrections, for samples diffused at dif- 
ferent temperatures. From these values the amount of diffusion 
occurring during heating and cooling was evaluated, assuming 
an exponential dependence D(T), which provided a second, 
more accurate value of D(T), This process was iterated 2 or 3 
times until it converged to constant values. Table I lists both 
the uncorrected and the corrected values of D. Error bars were 
obtained from the average of several EMPA scans, typically 
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Fig. 3- Extraordinary refractive index (open circles) and Zn dopant (trian- 
gles) profiles of an x-cut MgO:LtNbOj substrate after diffusion of a 200-nm 
ZnO film at 1000° C for lOh. The solid line is a complementary error function 
fitted to the experimental concentration profile. 



TABLE 1 

Zn Diffusion Coefficients For Different Temperatures 
And Substrates. The Corrected D Accounts For The 
Contribution Of Finite Diffusion During Warming And Cooling. 



Samples 


Temp. 


Time 


ZnO 

• thickness 


D(jim 2 fb) 


D (corrected, 


Mg0:LiNbO 3 
z-cut 


UOO^C 


1.4 h 


160 nm 


140±10 


1l9i8 




1100° C 


10 h 


160 nm 




24J±2 




950°C 


tQh 


1 fim 


S.7±2 


8S±2 




90Q°C 


100 h 


1 /im 


1.9±0JJ 


1.9±QJ 


WgOrLiNbO, 

X'CUl 


UOO^C 


1 h 


200 nm 


I38±I2 


110.8=1=10 


■ 


1000° C 


10 h 


200 nm 


27±6 


26.4±6 




900°C 


10 h 


200 nm 


2.4±0.2 


2,4=1:0.2 


TjNbQ 3 z-cul 


U00°C 


1 h 


1 ;<m 


2S0±26 


230±2I 




lOOO^C 


4h 


1 fim 




34±4 


• 


WC 


10 h 




9-6±1.4 


9.4±1.4 



2 to 4, taken at different positions on the same sample. They 
reflect some compositional nonuniformity across the samples 
as well as measurement uncertainties. Comparing the error bars 
of samples fabricated under similar conditions, it appeared that 
the dopant profile was more uniform for thinner ZnO films. 
This may be an important consideration in the preparation of 
waveguides with uniform index profiles. 

The samples used for D(T) measurements, which were 
fabricated with thick ZnO films, exhibit erfc profiles with a 
depth that varies as the square root of the diffusion time. This 
is an indication that, while second phase precipitation was 
present in these samples, the formation of this phase is clearly 
not a significant limiting step in the Zn transport process, at 
least within the accuracy of our measurements. Instead, the 
rate limiting step is diffusion of Zn in lithium niobate. The 
values of the diffusion coefficients D(T) measured from the 
diffusion of thick films should be applicable to the analysis of 
the diffusion of thin films necessary for waveguide fabrication. 

Fig, 4 shows the Arrhenius plot (ltiZ? versus l/T) for x- 
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Fig. 4. The Arrhenius (Jit D versus 1 jT) plot for the diffusion of ZnO. under 
conditions specified in the text, in z-cut and x-cut MgO:UNb0 3 between 900 
and 1]Q0°C 
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Fig. 5. Same as Fig, 4, for Zn diffused z-cut LiNoO^ between 950 and 
1 !00°C 

cut and z-cut MgO:LiNb0 3 . The dependence is Linear between 
900°C to 1100°C, as expected. The activation energy inferred 
from the slope of the curve is E 0 = 2.72± 0.2 eV. Within our 
measurement accuracy, we did not detect diffusion anisotropy 
in MgO:LtNb03, as shown by the results of Fig. 4. Fig. 5 
shows a similar curve for z-cut LiNb0 3 . The points falls on a 
straight line between the 950 and 1100°C, with an activation 
energy of Eq = 3.02 ± 0.3 eV. This measurement was not 
carried out for JC-cut LiNb0 3 . These high activation energies 
suggest that the diffusion mechanism of Zn in UNb0 3 favors 
substitutional over interstitial mechanisms, as the latter type is 
usually characterized by a smaller activation energy [20]. The 
diffusion of Zn is much faster than that of Ti, and fabrication 
times are correspondingly shorter. 
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IV. Index Characterization 
and Index-Dopant Relationship 

The index profiles, later used to establish the reJationship 
between the refractive-index change and the dopant concen- 
tration, were obtained by interference microscopy. For this 
purpose, the samples were sliced and polished as described 
in Section III. The samples, typically 1-mm thick, were 
placed in one of the arms of a Mach-Zehnder interference 
microscope operated with a 632.8 nm He-Ne laser, such that 
the light traveled through the length of the diffused region* 
The microscope produced a magnified image of the diffused 
region on which refractive index changes were mapped as 
distortions of the interference fringe pattern. Index profiles 
were inferred from the number of fringes of distortion at 
different positions along the depth of the diffused region. 
This method was generally preferred over prism coupling and 
IWKB analysis [21] because it is more easily applied to small 
samples with rough surfaces. In cases where the two methods 
were compared, the profiles agreed within 20%. 

Since the refractive- index change was typically small (less 
than 0.005), potential changes in the index from sources 
other than the Zn dopant had to be carefully evaluated. 
These sources include impurities in the ZnO film or in the 
furnace, departure from parallelism in the sample end faces, 
and Li in-or-out diffusion, which is known to change the 
extraordinary refractive index of lithium niobate [22], To 
make sure that no other process was affecting the index, a 
monitor sample of congruent LiNbOj (with no ZnO film) was 
routinely placed in the furnace and checked for contaminants 
by Jf-ray photoelectron spectrometry (XPS) and for index 
change by Mach-Zehnder interference microscopy. Within the 
resolution of both methods, we observed no detectable furnace 
contamination or Li in-or-out diffusion in the monitor sample. 
A slight wedge between the sample end faces and residual 
rounding of the edges can also be a source of error when 
interpreting the fringe pattern obtained from Mach-Zehnder 
interference microscopy. Michelson interference microscopy 
indicated that effects of rounding were negligible and that the 
contribution of any wedge to the measured refractive-index 
change was less than 1 x 10~ 4 . 

Studies of index profiles were carried out on samples dif- 
fused at 1100 and 1000° C for both x- and z-cut MgO;LiNb0 3 , 
and at 1000°C for z-cut LiNbOj. A typical index profile 
is shown in Fig. 3 for the extraordinary index of x-cut 
MgO:LiNb03. We estimate the error on the index measure- 
ment to be about S x 10~ 4 t and the spatial resolution to be 
2-3 fxm. The index change is maximum at the surface and 
equal to 3.7 x 10~ 3 . The half-maximum depth of the index 
profile is about 16 fim. The measured refractive-index profile 
resembles the Zn concentration profile measured in the same 
sample (see Fig. 3). Fig. 6 shows the dependence of the index 
change on the Zn concentration taken from the data of Fig. 3. 
The slope of a linear fit to the data is about 1.90 ±0.1 x 10~ 3 
per Zn mo\%. Fig. 7 shows a similar curve for another sample, 
also fabricated in x-cut MgO:LiNbO a but at 1100°C. For this 
particular sample, the index change appears to vary somewhat 
nonlinearly with Zn concentration. A linear fit to the data gives 
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Fig. 6. The extraordinary index change as a function of Zti concentration 
measured in an x-cul MgOrLiNbOj sample afisr diffusion of a 200 -nm ZnO 
film at 1000° C for lOh. 
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Fig. 7. Same as Fig, 6 t but for diffusion ai llOO'C for I h. 

a slope of about 1.77±0-12 x 10~ 3 per Zn mol%. More study 
needs to be done to define this dependence more precisely. 

Table II summarizes the index change per Zn mol% mea- 
sured on several samples representing a variety of materials, 
orientations and fabrication conditions. In both x- and z-cut 
MgO:LiNbC>3, the extraordinary index change is greater for 
diffusion at 1000°C than at HOCC. For a given material 
the index change appears to depend only slightly on the sub- 
strate orientation and diffusion temperature. These conclusions 
are based on a timited number of samples, and should be 
considered as tentative. 

For MgO:LiNb0 3 only the extraordinary index An e was 
studied, as the ordinary index change An 0 is generally too 
small to be measured accurately by interference microscopy, 
but was estimated to be on the order of 1/4 to 1/5 the value of 
An e . The ordinary modes were successfully excited by end- 
firing [11) but not by prism coupling, at least in *-cut samples, 
so that n 0 could not be measured by this method either. 
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TABLE U 

The Index Change Per Zs Concentration For 

DtFFE»ENT MaTCRJALS, Ow6NTATLONS AND TEMPERATURES. 



complementary error function (erfc) [23], which can be written 
as: 



Samples 


Temp. 


Time 


ZnO 
thickness 


Index 
change^ Zji 
mol%J 


I.MgO:LiNb0 3p x- 
cut 


HOO^C 


1 h 


200 nm 


1.77±0A2xlQ- 3 
(A*,) 


2, 


1000° c 


10 b 


* 


1.83±0.02xl0- 3 


3.MgO:LiNb0 3 ,x- 
cut 


1100°C 


1.4 h 


160 am 


1.22±0.04xl(r 3 


4. 


1000°C 


10 h 




1.90±(MQxlO- 3 


5. UNbO^-cut 


iooo°c 


4h 


200 nm 





For congruent LiNbCb the ordinary modes were easily 
excited by prism coupling. On the other hand the extraordinary 
modes were successfully excited by end firing but not by 
prism coupling. We also observed that the extraordinary index 
change is a nonmonotonic function which presents one or 
more extrema depending on the processing conditions. This 
behavior, which we believe may be related to local variations 
in lithium concentration, could be at the origin of the difficulty 
in prism coupling into this polarization. More studies are 
needed to elucidate this effect, which could be of interest to 
produce buried waveguides. This situation should be kept in 
mind when dealing with congruent LiNbOj. 

V. Comparing Calculated 
And Measured Index Profiles 

For the purpose of waveguide design, it would be useful to 
be able to predict the index profile of a waveguide from the 
fabrication conditions, i.e., from the film thickness, diffusion 
temperature, diffusion time, and the substrate material and ori- 
entation. To address this issue, we computed the concentration 
profile from diffusion theory, then multiplied this profile by the 
factor ATt/[Zn] (the measured ratio of index change An to 
zinc concentration [Zn] listed in Table II). We present results 
for Zn:LiNb03 in this section. The index profiles predicted are 
shown to agree within 30% error with the profiles measured 
from prism coupling and IWKB analysis. 

We considered the problem of one-dimensional diffusion 
of species A into solid B, with the assumption that (1) the 
solid solubility of A in B exists and is equal to C oy and 
(2) the diffusion coefficient of A into B is independent of 
concentration. We also assume that as long as the film is 
not depleted, the dopant concentration at the surface of solid 
B is pinned at the value of the solid solubility C a , This 
approximation for diffusion of thin ZnO films neglects the 
initial kinetics associated with dissolving Zn into LiNb03, 
but predicts results for waveguide samples with reasonable 
accuracy. Detailed modeling of the kinetic process is beyond 
the scope of this work. As discussed in Section II-B, in the 
present study this last assumption only applies when the ZnO 
film docs not exceed a prescribed thickness. In this case it 
is well known that the solution of the diffusion equation is a 



4* t t)«Co«rfc[-jL] (t<U) 



(1) 



where t is the total diffusion time and x Is the depth measured 
from the diffusion surface (i.e., the original interface between 
species A and B). This solution applies only until the film is 
completely depleted, which first occurs at time t = t\. The 
time ti is given by mass conservation, i.e., at t = ti the 
amount of material diffused into the solid per unit area is 
equal to the amount of material per unit area present in the 
film at time t = 0. Integration in x of (1) yields: 



.128C 0 

where Co is expressed in moie%, d is the initial film thickness, 
and k is given by: 

Pa 



PB M A 



(3) 



PA and pa are the density of species A and B, respectively, 
Ma M B are their respective molecular weight, and N is 
the number of tons per molecule. For the case under study 
A/=l and k = 2,23, 

For t > ti, the solution of the diffusion equation is known 



(4) 



where t2 = t — ti- In general, this solution can not be cast 
in a simple closed form expression and the integral has to be 
calculated numerically. As is well known, for t% > t\ this 
solution approaches a Gaussian given by: ■ 

Kd 



<T 4^1,4^ (f 2 > ti) (5) 



The 1/e width in (5) is w G = 4L> V /(£ 1 + ^ exact 
solution (4) converges toward the approximate solution fairly 
rapidly. For t 2 = 1 1 the exact width differs from the approx- 
imate width wo by only 11%, and for t% = 4£i by less than 
4%, In practical diffusion runs, ±2 is usually longer than t\ t 
so that (5) provides a convenient and accurate approximation 
of the dopant profile. 

The solid solubilities of Zn in LiNb0 3 and MgOiLLNbOj 
are necessary inputs for this calculation. In an ideal system, 
we would expect the solid solubility to depend on tempera- 
ture, but at a given temperature the Zn concentration at the 
surface should be pinned at that value as long as the film 
is undepleted. In the ZmLiNbOs system, in which the ZnO is 
not in thermodynamic equilibrium with the substrate materials, 
we obtained an effective solid solubility by measuring the Zn 
surface concentration with EMPA in several waveguides (with 
no detectable second phase precipitation, i.e., starting from 
ZnO film thicknesses smaller than the bounds specified in 
Section II-B). For a given temperature and substrate material, 
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the highest measured value provided a lower bound value of 
the solid solubility. For example, in a waveguide fabricated 
at 1100*C in x-cut MgO:LiNb0 3 the highest Zn surface 
concentration was equal to 3.5 ± 0.6 mo!e%. We assume C 0 
lakes this value for all temperatures between 1000 to 1100°C, 
as no data was available for 1000°C. For LiNb0 3( the highest 
surface concentration observed and the value of C a used was 
5.4 ±0,5 mole%. These two values were consistent with mass 
conservation, assuming that the film was just depleted at the 
end of these runs. 

The diffusion temperature profile 7X0 is not rectangular in 
practice, as assumed in the model. Instead, the furnace is grad- 
ually heated to the diffusion temperature Tq, the temperature is 
then maintained at T 0 for a duration t t and finally the furnace is 
gradually cooled down. As mentioned earlier, a finite amount 
of diffusion takes place during heating and cooling cycles. The 
real cycle can be replaced by an equivalent rectangular cycle 
(i.e., with infinitely fast heating and cooling) producing the 
same amount of diffusion, provided time t is replaced by an 
equivalent diffusion time f efl which can be shown to be given 

where k is the Boltzman constant, E 0 is the activation energy 
of the dopant diffusion in the substrate, and 0+ and /?_ are the 
healing and cooling rates, respectively. The correction term in 
(6) becomes negligible when t is large. For the heating and 
cooling rates used in this study, in the case of ZnO diffusion 
in lithium niobate the correction is negligible when t exceeds 
about 10 h. 

To verify the applicability of this model to Zn-diffused 
waveguides, we compared the measured index profiles of six 
waveguides with profiles calculated as follows: 1) calculate the 
corrected diffusion coefficient D at the diffusion temperature 
from Table I, 2) calculate both the effective diffusion time 
k// (6) and t { (2), 3) calculate the dopant profile, using (1) 
if t vff < h or (4) if t eff > ii (or (5) if t e// > 2*0, and 4) 
multiply Ihe dopant profile by the An/[Zn] factor (Table IT) 
to obtain the index profile. 

The waveguide mode effective indices and mode turning 
points were measured by prism coupling and IWKB analysis. 
In the IWKB analysis a curve fitting routine was used to fit 
the experimental index profile to either a erfc or a Gaussian 
(depending on whether t was smaller or larger than ti) and 
determine the index at the surface. Error in the effective 
index measurement was ±0.0003. For a typical three mode 
waveguide fabricated in this research, the uncertainty in the 
assignment of the surface index in IWKB analysis calculation 
and the mode turning points are about ±0.0006 and ±0,6 m, 
respectively. 

Fig. 8 shows the measured index profile and the calculated 
index profiles for the extraordinary polarization of a waveguide 
in MgO:LiNb0 3 . The fabrication conditions were SO-nm ZnO 
film diffused at 1000°C for 1 h in a *-cut sample. Including the 
finite ramping time, the effective diffusion time was calculated 
to be 72 min. The calculated value of *i (2) was 45 min. In 
this case, t 2 = 72-45 = 27 min. The solid and dashed lines in 
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Fist 8. Extraordinary itvdex data rawsured by IWKB analysis on a waveguide 
fabricated in x-cui MgO:LiNbO, by diffusion of a 80-n/n ZnO film at 100ff> C 
for 1 h- The solid curve is ihe exact calculated profile (4) and the doited line 
the approximate calculated profile (5). 




FIk 9. Similar as Fig. 8, for z-cut LiNbOj (ordinary index data) after 
diffusion of a 100-nm ZnO film at 1000°C for 45 min. 



Fig. 8 represent the profiles calculated according to the exact 
solution (4) and the approximate solution (5), respectively. 
The agreement between the exact profile and the measured 
profile is reasonable. The measured index value at the surface 
is somewhat higher than predicted, which could be due to 
uncertainties in the measured value of AnfiZn], as discussed 
in Section IV. 

Similar results are presented in Fig. 9 for a waveguide in 
LiNb0 3 (100-nm thick ZnO film, r-cut sample, 1000°C, 45 
min). Including the finite ramping time, this corresponds to 
an effective diffusion time of 56.6 min. The calculated value 
of ti (2) was 23 min. In this case, t 2 = 33 min, the ratio 
t 2 /h = 1.4 is large enough that the exact profile resembles 
a Gaussian. Again the theoretical profiles are in reasonable 
agreement with the measured index profile. 
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Fig. 10. Same as Fig. 9, with a diffusion time of 12 min. for two waveguides. 
Both waveguides were fabricated under the 4am e conditions but in different 
runs- The solid curve is the calculated index profile calculated from (1). 

Fig. 10 shows similar results for two waveguides as in 
Fig, 9 except that the diffusion time was decreased to 12 
min. These two waveguides were fabricated under the same 
conditions but in different runs. Including the finite ramping 
time, the effective diffusion time was calculated to be 23.6 
min. In this case ti — 23 min. so that the film was just 
depleted at the end of the diffusion process and (1) was used to 
generate the theoretical profile of Fig. 10. The surface index 
for the 1WKB fit was determined by fitting to an erfc. For 
one of the waveguides the measured and calculated profiles 
agreed very well, to better than a few percent. For the second 
waveguide the agreement was fair, the actual profile depth was 
approximately 20% larger than the calculated depth. 

The difference between the two data sets obtained from two 
different waveguides illustrates the limits of the reproducibility 
of our fabrication process. We believe this sample-to-sample 
variation is largely due to variations in film thickness. These 
variations change the total amount of Zn diffused into the 
sample, and therefore the total area under the index profile (see 
Fig. 10). Also, we speculate that film thickness variations may 
influence the total rate of Zn transport through the following 
process. While diffusion rate in single phase UN0O3 is similar 
for both thin and thick ZnO film samples, the initial kinetics 
associated with Zn dissolution into L1N0O3 might require a 
period of time f sp . In this ZnO-LiNb03 nonequilibrium system, 
may differ for films of different thickness. As discussed 
earlier, for typical diffusion times (30-60 mirr) this initial 
transient is not a significant rate limiting step and has a 
negligible effect on the final waveguide depth. However, in 
the case of the two waveguides of Fig. 10 the diffusion time 
was unusually short (t — 12 min), and the initial transient may 
influence the final waveguide depths. If depends on the film 
thickness, in this regime of very short diffusion times sample- 
to-sample variations in film thickness may lead to variations 
in waveguide depth (see Fig, 10) 

The index profile is fairly sensitive to the fabrication con- 
ditions, so that the accuracy of our model in predicting a 



given waveguide profile is limited by the accuracy of the input 
parameters. For example, in the case of Fig. 9, a nearly perfect 
fit to the experimental profile is obtained when increasing the 
ZnO film thickness by 15% and decreasing the value of D by 
15%, These variations are within the experimental errors in 
our measurements of these two quantities. The predictions of 
the model presented here should become more accurate when 
the process is better characterized. 

Similar agreement was obtained for three other waveguides. 
Together with the examples illustrated in Figs. 8-10, these 
comparisons indicate that our model provides reasonable pre- 
dictions of the index profile, generally within 40% for the 
surface index change and 20% for the 1/e depth. More detailed 
measurements are needed to refine the values of the model 
parameters, i.e., the diffusion coefficients and the An/[Zn] 
coefficients. 

VI. Conclusions 

ZnO diffusion in MgOzLiNbOj and congruent LiNbOj, and 
the index profile of waveguides fabricated by this process have 
been characterized. Earlier reports established the usefulness 
of this new type of waveguide, which guides both polarizations 
with low propagation losses and, in the case of MgOiLiNbOs, 
exhibits high resistance to photorefractive damage. In the 
present study it was found that using ZnO instead of metallic 
zinc as a source and maintaining the ZnO film thickness 
below a prescribed value were important steps in avoiding 
second phase precipitation and maintaining the surface quality 
required for low loss waveguides. The refractive-index change 
produced by the Zn dopant varies approximately linearly 
with Zn concentration, with only a slight dependence on 
substrate orientation and diffusion temperature. Simple ex- 
pressions from linear diffusion theory predict the waveguide 
index profile from the fabrication conditions. Comparison 
between calculated and measured profiles for a variety of 
waveguides showed that these model predictions are fairly 
accurate. Further studies of the dependence of the index 
change on Zn concentration and of the extraordinary index 
change in congruent LiNbOj are needed to improve the 
understanding of these waveguides. Diffusion studies using 
suitable dopant sources within the tie triangles containing 
LiNbOj can potentially provide useful information about the 
Zn transport. ZnO diffusion in lithium niobate appears to be 
an interesting process which has the potential to extend the 
spectrum of available waveguide devices. 
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* „r t« hnc h^n demonstrated on a -'/-cut substrate at 0.6328-;tm wave- 
Th mmsured results show that the »a v e»uide s s u p p or t conhned waveguide at the n e a r -IR wav e l e ngth region due- 



reported. The measured results show ihai ihe W a>«= S u.Jts support 
either a single extraordinary polarization or both extraordinary 
and ordinarv polarizations depending on the fabrication process 
parameters/For the single extraordinary polarization waveguide, 
ihe measured propagation loss at L32- M m wavelength is 0 9 dB/cm 
and the best measured polarization extinction ratio is 44 dB at a 
distance of L5 cm from the input end, which are quite good for 
being a waveguide polarizer. Moreover, the voltage-length product 
measured for the ZI Mach-Zehnder modulator shows that the sub- 
strate electrooptic coefficient is not degraded. 

Index Terms— Lithium niobate, polarization-dependent loss 
measurement, polarizer, Zn indifTusion* 

I. introduction 

SINGLE polarization waveguides are widely used for the 
fabrication of integrated optical sensors because a single 
polarization wave is essential to increase the optical sensitivity 
[1]. Moreover, a combination of a single and both polariza- 
tion waveguides can be a polarization splitter with large split- 
ting extinction ratio and fabrication tolerance [2]-[4]. To date, 
metal-clad Ti indiffusion (TI) and proton-exchange (PE) are two 
well-known methods for making single polarization waveguides 
on LiNb0 3 [5], [6]. In the metal-clad waveguide, the TM-wavc 
propagation attenuation, depending on the buffer layer thickness 
and chosen metal-clad, is much larger than that of the TE~wave, 
After propagating a certain distance, only the TE-wave is left 
in the waveguide and it will be slightly affected by the cladding 
lavers. As to the PE waveguide, the change in the extraordinary 
index in e is positive while the change in the ordinary index 
An. is negative. Thus, it supports only an extraordinary po- 
larization wave. An annealing process is often used to reduce 
the propagation loss and restore the electrooptic coefficients. 
Therefore? good-performance devices have been successfully 
demonstrated on the r- and x-cut substrates [6], [7]. However, 
proton-exchanged waveguides suffer from surface damages and 
the guiding properties are consequently degraded on y-cux sub- 
straps [S]T This may cause serious problems in surface-wave 
acoustic-optical applications, since surface acoustic wave de- 
vices frequently require //-cut substrates. Although a tunable PE 
directional coupler with a Si0 2 cladding and thermal annealing 
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coniined waveguide at ihe n e a r -IR wav e l e ngth region duc-xc 
the PE depth and concentration limitations for avoiding surface 
damage. 

Recently, Zn indiffusion (ZI) waveguides have been studiec 
extensively to make optical devices in LiNb0 3 due to les: 
susceptibility to photorefractive damage in comparison t< 
TI waveguides [10}. Also, the Zn atoms have larger diffu 
sion coefficients than those of the Ti atoms. Therefore, Zi 
diffusion needs no additional out-diffusion suppression an* 
saves the time for diffusion, which provides flexibility for th 
fabrication of optical waveguides. The guiding characteristic 
at 0.6328-/im wavelength of a planar waveguide have bee 
discussed [10], [11]. However, the variations of guided mod 
profiles and the polarization-dependent losses of a chanm 
waveguide fabricated with different process parameters ha\ 
neverbeen addressed at the L32-/jm wavelength. 

In this paper, the polarization-dependent loss rneasurernen 
of ZI waveguides on a y-cvu ^-propagation LiNb0 3 substra 
are demonstrated. The results show that the ZI wavegmc 
supports a single extraordinary or both extraordinary and o 
dinarv polarizations depending on process parameters. In sue 
sinsJe extraordinary polarization waveguide, the propagatic 
loss is 0.9 dB/cm and the best polarization extinction rat 
measured at a distance of about 1.5 cm away from the inp 
end of the waveguide is 44 dB. Therefore, it is quite suitat 
for using as a waveguide polarizer. Moreover, the fabricat 
ZI Mach-Zehnder modulator shows that it has a nondegrad 
electrooptic coefficient r 33 . On the y-cut substrate, a ZI wai 
suide polarizer has the advantage of better performance a 
simpler fabrication process than that of a PE or a metal-clad 
waveguide polarizer 



IL Experiments 

First, a Zn thin film is deposited on the substrate by then 
evaporation. As Zn atoms are difficult to be leeched directly 
the surface of LiNbO^ substrate, a 50- A Ni film is predeposi 
onto the LiNbOj substrate in order to increase the adhesion 
channel waveguide and a Mach-Zehnder modulator patterns 
then formed by the lift- off technique. The waveguide widtl 
S u.m, the r -branch angle is 2°, and the gap between two ai 
of the Mach-Zehnder modulaior is 16 |im. The samples v. 
placed on a covered alumina crucible, which was then place* 
a high-temperature oven. The thermal cycle, which is 
in dry air, typically involves heating at a rate of about 10 L/i 
up to the diffusion temperature, maintaining at that tempera 
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for a certain lime, and cooling at a rale of 10 p C/min down to 
the room temperature. 

After the thermal diffusion, the substrate end faces are cut 
and polished 10 allow end butt coupling. A schematic diagram 
of the polarization-dependent loss measurement setup is shown 
in Fig, 1. A He-Ne laser of wavelength 0.6328 //m is used for 
optica! alignment. The polarization of an incident Nd:YAG laser 
light of wavelength 1 .32 //m is controlled by rotating the polar- 
izer. The angle is varied from 0 = to 360' relative to the ;i/-axis 
with a step of 10°, The incident light is in the TM polarization 
when the angle is set at 0° or 3 SO*, and in the TE polarization, 
when the angle is set at 90 c or 270°. For comparison, an attenu- 
ator is used to keep the incident power constant at different po- 
larizations. The light is divided into a reference and a measure- 
ment light after passing through the beam splitter. The measure- 
ment light is coupled into the front end face of the waveguide 
w ith a 40x lens and the output beam is imaged onto an InGaAs 
phoiodeiector (or a charge-coupled device camera) also with a 
40x lens. Both reference and measurement lights are detected 
by a dual-channel power meter at the same time. The measured 
output mode profile is displayed on a video monitor. 

For the fabrication of a Mach-Zehnder modulator, a pair of 
aluminum electrodes of thickness 3000 A, length 6 mm, and 
gap 14 //m is formed forelectrooptic modulation. The switching 
voltage measurement setup is also shown in Fig, 1. A triangular 
voltage signal from the function generator is applied across the 
electrodes. The input and measured signals were both shown 
on an oscilloscope which were further used to calculate the 
switching voltage. 

III. results and Discussion 

The mode profiles are dependent on the strip film thickness 
t\ waveguide width IT', diffusion temperature T, diffusion time 
/, and operation wavelength A. As single-mode waveguides are 
essential for most device applications, only the single-mode 
waveguides are further discussed latter on. The TE (extra- 
ordinary polarization) and TM ^ordinary polarization) mode 
profiles versus different process parameters are shown in Fig. 
2. Consider case .4 as illustrated in Fig, 2. The TE and TM 
modes show that they have similar guiding profiles. In case 
B. the TE mode size is slightly larger than that of the TM 
mode. In case C y the TM mode is not well-confined in the 
waveguide and most of the power is radiated into the substrate. 
However, only the TE mode is guided in case D- The above 
experimental results show that the Zl waveguide supports 
a single extraordinary or both extraordinary' and ordinary 
polarizations depending on fabrication process parameters. 
Besides, the single extraordinary potari7ation waveguide is 
easier achieved than the single ordinary polarization waveguide 
due to Ait t > ^tt.,, particularly when fabricated at a lower 
diffusion temperature [11]. Also, well-confined TE and TM 
waveguides are more difficult to be obtained when the diffusion 
temperature is below 900 C. 

In order to further siud\ the polarization-dependent loss, the 
normalized output power versus diffcreni polarizations with 
some process parameters is shown in Fig. 3. As the incident 
power is kept constant for any polarizations, the polarization 
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Fig. I. A schematic diagram of (he polarizaiion-dependem loss and switching 
voltage measurement .setup. 
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extinction ratio A is simply defined as = 101o£( Ptr/Pjw K 
where /Ve -^tm are the ouipui powers of TE and TM 
modes, respectively. Jn case the best 7 obtained is 44 dB at a 
distance of about 1.5 cm away from the input end of the wave- 
guide and the propagation loss measured by a cut-back method 
is 0.9 dB/cm. Fig. A gives the modulation characteristics curve 
of a Mach-Zehnder modulator with its waveguides fabricated 
as those of case Z~>. The measured switching voltage is 17.5 V. 
the corresponding voltage- length product is 10.5 V-crn. the 
calculated overlap integral factor by the reported method |I2] 
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1.32-/*m wavelength fabricated with different process param- 
eters is firstly and successfully reported. The results show 
that the ZI waveguide supports a single extraordinary or both 
extraordinary and ordinary polarizations depending on the 
process parameters. For the single extraordinary polarizarion 
waveguide, the best 7 measured is 44 dB and the propagation 
loss Is 0.9 dB/cm. These are all good enough for making 
a waveguide polarizer. The fabricated ZI Mach-Zehnder 
modulator also shows that it has a nondegraded elecirooptL 
coefficient 7-33. With the excellent properties of process-de- 
pendent polarization, more flexibility on waveguide device 



fabrication is of great interest in the near future. 
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Fig. 3. The normalized output power versus polarization for waveguides 
fabricated with different process parameters. 




Fig. 4. Measured sinusoidal output signal (2 nVWdiv) of a zinc indiffusion 
Mach-Zehnder interferometer with an input saw-tooth wave (5 V/div). 



is 0.45. From the above data, the derived r 33 coefficient is 
about 30.8 pm/V, which shows that the coefficient of a ZI 
waveguide is not degraded under the specific process condition 
presented in this work. 

IV. Conclusion 

A polarization-dependent loss measurement of ZI waveg- 
uides on a y-cut, .r-propagation LiNbOs substrate for the 
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Fabrication of Zn:LiNb0 3 Waveguides by Diffusing ZnO in Low-Pressure Atmosphere 
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Fabrication of ZmLiNbCK planer and channel waveguides by diffusing ZnO in low-pressure air is presented. The process 
is simple and does not require special apparatus. Pressure control is crucial to ensure modest indiffusion with suppressed out- 
diffusion. The waveguides support both ordinary and extraordinary polarization modes, and arc expected to be phoiorcfractive- 
damage resistant. The diffusion temperature considerably lower than that for Ti-tndiiTusion is advantageous for applications to 
quasi-phase matched nonlinear optic devices. Single-mode guides of — 1 dB/cm propagation loss for 1.5 jjm wavelength band 
were obtained. 

- KEYWORDS: ■ opt i c al wav e guid e , l i thium n i ob a te , quas i p h a se m a t c hing ,. non li ne a r optic s. . d i ffu s ed . waveguide 



1. Introduction 

There has been much research and development inter- 
est in quasi-phase-matched (QPM) nonlinear-optic (NLO) 
devices and rare-earth doped lasers in LiNb0 3 wave- 
guides. Annealed-proton-exchanged (APE) waveguides' - y 
are widely used for such applications. The APE waveguides, 
however* have a drawback that they support only extraordi- 
nary guided modes, and it may pose limitations in implemen- 
tation of polarization-independent NLO wavelength conven- 
ers and integration of NLO t laser and electrooptic devices, Ti- 
indiffused guides, 3 * which support both ordinary and extraor- 
dinary modes, involve the photorefractive damage problem. 4 * 
Fabrication of periodically-domain-inverted grating for QPM 
in Ti:LiNb0 3 ^uide requires a complicated process 5) to solve 
the problem' caused by unwanted domain inversion during the 
high-temperature Ti diffusion process. 

An alternative waveguide is Zn-indiffused waveguides, 
which support both polarization modes and are considered to 
be photorefractive-damage resistant. 6- ^ Planar guides were 
fabricated by diffusing ZnO ai 1000-1 100°C into z-cut un- 
doped and MgO-doped LiNb0 3 crystals in Li-rich atmo- 
sphere, and were characterized at 0.63 tim wavelength. 6 * 
Er-doped lasers have been demonstrated using channel guides 
fabricated by diffusing ZnO into z-cut Er:MgO:LiNb0 3 at 
930°C for 0.5 h in air, although the detail of the atmosphere 
was not presented. 7 * Another feature of Zn:LiNb0 3 guides 
is that they can possibly be fabricated by diffusion at tem- 
peratures lower than that for Ti- indiffusion. In fact t planar 
guides were fabricated by metallic Zn diffusion from vapor- 
phase at 600-700°C in low-pressure nitrogen atmosphere us- 
ing a specially designed apparatus, and were characterized 
at 0.44-0.63 Atm wavelengths. 85 Channel guides for 1,3 j^m 
wavelength were fabricated by diffusing Zn film into y-cut 
LiNb0 3 . 9) However, techniques for fabrication of Zn:LiNb0 3 
guides have not been fully established. 

In this letter, we present fabrication of Zn-indiffused planer 
and channel waveguides for up to 1 .5 Mm wavelength band in 
z-cut LiNb0 3 by diffusing ZnO in low-pressure atmosphere. 

2. Planar Waveguides 

As a preliminary experiment, diffusion of ZnO was tried 
with various atmosphere. Congruent optical-grade z-cut 
LiNb0 3 crystals were used, and thin film of ZnO of 280 nm 
thickness was deposited on the -f-z surface by RF magnetron 
sputtering in Ar/02 (1:1) atmosphere using a ZnO target. An 



ordinary quartz-tube electric furnace was used for heat treat- 
ment at 930°C for 0.5-24 h. The surfaces of crystals heated in 
(dry and wet) 0 2 and air at 1 atm were rough and were covered 
with nonuniform residual which seemed to have been yielded 
by ZnO-LiNb0 3 reaction, and therefore waveguide was not 
obtained with such samples. This result is in contrast with 
the successful waveguide fabrication in Er:MgO:LiNb0 3 by 
similar process. 7 ^ The surface of crystals heated in quartz tube 
evacuated by a rotary pump, on the other hand, was clear, and 
planar guides were obtained. Waveguides obtained by diffus- 
ing 175 nm thick ZnO at 930 6 C for 0.5 h, for example, sup- 
ported transverse electric (TE) and transverse magnetic (TM) 
modes at 0.46, 0.63, 1.06 and 1.54 tim wavelengths. The 
refractive-index increment at the surface estimated from on- 
line measurement at 0.46 /um wavelength was 0.006 for both 
ordinary and extraordinary waves. Almost same results were 
obtained using an lR-heating rapid thermal annealing (RTA) 
furnace. 

3. Channel Waveguides 

To examine further the dependence of the diffusion on 
the process parameters and find conditions for fabrication of 
waveguides for device applications, channel guides were fab- 
ricated. The ZnO film of 220 nm thickness was patterned 
into 7 Atm wide stripes by lift-off technique. The air pres- 
sure in the quartz tube was monitored with a Pirani gauge and 
was controlled by a throttle valve. The crystals were slightly 
darkened due to reducing by heating in oxygen-poor atmo- 
sphere, depending upon the temperature, time and pressure. 
The crystals heated at lower pressure were darker. The origi- 
nal transparency was recovered by annealing in O2 0 atm) at 
650°C for 5 h. Examination of the crystal surface after etching 
in HF-HNO3 showed no evidence of unwanted ferroelectric- 
domain inversion. The both ends of the waveguides were pol- 
ished for end-fire coupling. The waveguide length was ap- 
proximately 10 nam. The waveguiding characteristics were 
examined at a few wavelengths. The results are summarized 
in Table I. 

It was found that the diffusion depends significantly on the 
pressure. For guides diffused at 930°C in pressures at or lower 
than ITorr, channel guiding was weak or not observed, and 
evidence of planar guiding of TM mode only was observed 
as a result of increase of extraordinary index due to Li 2 0 
outdiffusion enhanced in low pressure. For guides diffused 
at lOTorr, on the other hand, only channel guiding was ob- 
served. The fact that no channel guiding was observed for 
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Table 1. Result of ihe optical waveguiding test of the fabricated channel guides. 


Sample 


Diffusion 


Evalution wavelength, polarization 


Tcmprature, Time. Pressure 


0.63MmTM/re 


1.54/imTMTE 


#1 


930°C 0.511, lTorr 


p/n 


p/n 


#2 


930°C 0.5 h, lOTon- 


s/s 


s/n 


#3 


S90°C 0.5 h t 0.6Torr 


m/m 


s/s 


#4 


890 & C. 0.5 b, lOTorr 


m/m 


s/s 


#5 


890*C. I.Oh, lOTorr 


m/m 


m/s 


#6 


S90°C, 2.0b. lOTorr 


m/m 


m/s 



p: planar guiding s: single- mode guiding m: muUimode guiding n: no suidine. 




Fig. 1 . Near field patterns of the guided modes at 1 ,54 /im wavelength ot> 
serbed in channnel guide #4. (a) TM-Hke mode, (b) TE-Uke mode. 



samples diffused at higher temperature and in lower pressure 
may be explained by insufficient index increment due to fast 
(too deep) diffusion , although further experimental study is 
required to clarify the reason. These results indicate that for 
fabrication of device-quality cannel waveguides pressure con- 



trol is crucial to ensure modest indiffusion with suppression 
of outdiffusion, support the observation that the diffusion of 
Zn is enhanced by oxygen vacancies,*) and also suggest the 
possibility of fabrication at even lower temperature. 

The guides #3 and #4 supported single mode at 1.54 fim 
wavelength. The mode confinement was stronger in #4 than 
in #3. Figure 1 shows the mode patterns of guide #4 dif- 
fused at 930°C and lOTorr. The full width at half-maximum 
(FWHM) mode size was 4.3 fim (depth) x 7.0 jxm (width), 
and 7.7 Atm x 14 A<m for TM-like and TE-like modes, respec- 
tively. The propagation loss was measured as approximately 
1 dB/cm for both modes by the waveguide Fabry-Perot inter- 
ference method using a temperature-tuned distributed feed- 
back (DFB) laser 

4. Conclusion 

We demonstrated successful fabrication of Zn-indiffused 
LiNbCh waveguides by diffusing ZnO in low-pressure air. 
Control of the pressure is crucial for obtaining device-quality 
guides with good reproducibility. The process is simple and 
does not require a special apparatus. The guides support both 
polarization modes, and are expected to be photorefractive- 
damage resistant. Further experimental study is required to 
optimize the diffusion conditions and clarify the photorefrac- 
tive characteristics. The diffusion temperature considerably 
lower than that for Ti-indiffusion is advantageous for appli- 
cation to QPM-NLO devices. Most recently, we actually fab- 
ricated QPM second harmonic generation devices and differ- 
ence frequency generation devices, and obtained encouraging 
results. The results will be published elsewhere. 
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LiNb0 3 optical waveguides by Zn diffusion from vapor phase 
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Fabrication of optical waveguides on LiNb0 3 by Zn diffusion from the vapor phase has been 
demonstrated, using a specially designed diffusion apparatus. The planar waveguides on z-cut 
samples support both TE and TM propagations, indicating an increase of refractive indices for n Q 
and n t , Energy dispersive x-ray analysis measurements confirm the diffusion of Zn into the crystal, 
and analysis of the content of Zn as a function of diffusion depth for samples treated at different 
temperatures indicates much higher diffusion coefficients than those obtained for Ti diffusion in 
LiNb0 3 . The technique is versatile, easy to operate, and allows fabrication of waveguides at 
temperatures as low as 600-700 °C, thus preventing Li outdiffusion. O 1995 American Institute of 
Physics. 



LiNb0 3 is an important inorganic key material for opto- 
electronics, which combines high optical nonlinearities with 
the capability of incorporating high concentrations of laser 
active ions. Several methods are now widely used for fabri- 
cating waveguides in this materia!, 1 " 3 the most popular being 
Ti diffusion, proton exchange, and ion implantation. Each 
one of the developed techniques has some limits of applica- 
bility: proton exchange only increases in the extraordinary 
refractive index; ion implantation requires the use of ion ac- 
celerators. Even though, Ti diffusion does not have these 
drawbacks but the process is performed at temperatures rang- 
ing from 950 to 1200 °C, 4 and at such high temperatures the 
Li-outdiffusion problem (which could form undesirable sur- 
face waveguides competing with channel waveguides 5 ), must 
be addressed. Several methods for suppressing outdiffusion 
have been proposed, 6 which add complexity to the technique. 

To overcome such difficulties, low temperature diffusion 
is desirable; this requires the use of elements that have high 
diffusivity and low activation energy, and that, upon diffu- 
sion, could cause an index increase and do not introduce 
excessive absorption at the operating wavelength. 

Under these requirements, metallic Zn appears to be a 
potential candidate. In fact, Zn diffusion from the vapor 
phase has been demonstrated to produce low loss optical 
waveguides in LiTa0 3 . 7 The procedure consists of loading a 
bare, cleaned substrate together with Zn wire or granules in a 
quartz ampoule, which is evacuated with a diffusion pump 
and sealed. The diffusion is then performed by heating the 
ampoule in a tube furnace. Nevertheless these authors re- 
ported no success in waveguide formation in LiNb0 3 due to 
severe pitting and surface damage on the substrates. 8 

In the present work we report a successful method to 
produce monomode and multimode optical waveguides in 
LiNb0 3 by metal Zn diffusion in a controlled way at tem- 
peratures as low as 600 °C. This has been achieved by dif- 
fusing the metal Zn from the vapor phase in a specially de- 
signed diffusion apparatus. We also report measurements by 
energy dispersive x-ray analysis (ED AX) of the dopant dis- 
tribution in Zn-diffused LiNbOj waveguides, as well as the 
optical characterization of their index profiles. These data are 
used to obtain values of the diffusion coefficient of Zn in 
LiNbO-, their temperature dependence and the relationship 
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between the refractive index change and the Zn concentra- 
tion. Comparison between diffusion coefficients of Zn, ZnO, 
and Ti establishes the potential advantage of this method for 
fabricating waveguides in LiNb0 3 . 

Single crystals of congruent LiNbOj were grown in our 
laboratory by the Czochralski method from grade 1 Johnson- 
Mathey powder. The initial preparation of the samples in- 
volved cutting slices from the boule and polishing the crys- 
tals to optical quality. The finished samples were all of the 
same cut (z cut) and had final dimensions of 2X6X8 mm 3 . 
The diffusion apparatus used in this work is based in a Van 
Doom heatpipe, traditionally used in the additive coloration 
of alkali halides. The setup is described in the literature, 9 and 
consists basically of a stainless steel tube previously cleaned 
at high temperature, where the metal is deposited on the 
bottom, and the sample is placed on a basket, Stablizied 
347-type stainless steel was used to prevent corrosion in the 
bottom (welded zone) of the tube. Water cooling in the upper 
part and a conic condenser were used to confine the metal 
vapor. The upper part of the tube is a standard NW-20 KP 
flange connection through which the tube is either evacuated 
or filled with nitrogen at the desired pressure, also providing 
space for the thermocouple feedthrough and the crystal con- 
tainer fixing. The tube is finally introduced in an oven, the 
heating and cooling rates being 10 °C/min, to avoid thermal 
shock to the crystals. 

Zn diffusion profiles were determined by EDAX mea- • 
surements. Prior to EDAX characterization, Zn-diffused/ 
samples were first sandwiched between two glass slides, and * 
the end perpendicular to the diffusion surface was polished. 
Glass slides were used to prevent rounding of the edge of the 
diffused region during polishing, which would otherwise up- 
set EDAX measurements. Optical characterization of the 
samples includes determination of the refractive index pro- 
files of the planar waveguides, computed by fitting the mode 
spectra measured at \=442, 488, 514, and 633 nm with ap- 
propriate Gaussian index profiles. 

Z-cut samples of LiNb0 3 were treated at 500, 600, and 
700 °C during 0.5, 4, and 4 h respectively (hereafter labeled 
as Nos. 1, 2, and 3), using a nitrogen buffer pressure of 10 
Torr, with ~2 g of metallic Zn. This pressure is adequate in 
this range of temperature to provide high enough Zn vapor 
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TABLE I. Number of observed dark modes for sample Nos. 2 and 3 at 
different wavelengths. 
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FIG. 1. Zn cuncenlritiun in atomic percentage, measured by EDAX tech- 
nique, as a function of depth normal to ihe diffused surface for LiNbOj 
crystals treated at different temperatures in Zn metallic vapor atmosphere 
(lines are only a guide for the eyes), 

surrounding the sample, thus allowing the diffusion. After 
the diffusion, the samples showed a dark gray color coming 
from color center formation during this process. Subsequent 
annealing in air removed these color centers and the samples 
became transparent. 

We have followed the metal diffusion through the crys- 
tals by means of the EDAX technique. After diffusion, one 
end of the samples was repolished, and EDAX scannings of 
the edges were performed. Figure 1 shows the Za concentra- 
tion (at. %) plotted as a function of diffusion depth for the 
three samples. As can be observed, sample No. 1 has a peak 
value in %Zn of *- 1 .2%, with a penetration depth (computed 
at l/e 2 of its maximum value) of —0.2 /xm. Samples Nos. 2 
and 3 have similar values for a Zn concentration near the 
surface of --4.7%, and diffusion depths of —0.9 and —3,0 
fim, respectively. This concentration at the surface for 
sample Nos. 2 and 3 is close to the highest surface concen- 
tration observed in LiNb0 3 :Zn waveguides formed from 
ZnO at higher temperature 10 (5.4%), indicating similar val- 
ues of the effective solid solubility for both fabrication pro- 
cedures. 

The calculated values for the diffusion coefficients are 
£>(600 °C)~*0.06 /xm 2 /h and £(700 °C)~0.6 jLtmVh. These 
coefficients are about two orders of magnitude higher than 
the diffusion coefficients of Zn from ZnO source at these 
temperatures, calculated by extrapolating data given by 
Young etal., 10 using Z>,> = 4,2Xl0 n ^m 2 /h, E Q 
= 3.0 eV. Although the origin of this discrepancy is uncer- 
tain, a possible origin could be connected with the rate of 
oxygen vacancies annhilation depending on the atmosphere, 
as suggested by Burns 11 for Ti indiffusion in LiNb0 3 . More- 
over, the value of D at 700 °C is comparable to the reported 
diffusion coefficient for Ti at 1100 °C, U indicating that at a 
given temperature, the diffusion of Zn from the vapor phase 
is much faster than that of Ti. This behavior is similar to that 
found in Zn vapor diffused Zn into LiTa0 3 , where the 
authors 8 report a value of D{800 °C)=0.7 /xm a /h, which may 
be compared with £>(1 200 6 C)= 0.29 /xm 2 /h for Ti 
diffusion, 11 

Optical waveguiding properties of the samples were then 
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examined by dark mode experiments after removing the 
color centers. Sample No. I, treated at 500 °C, did not show 
any dark mode either for TE or TM polarizations. The fact 
that this sample did not reveal optical waveguiding behavior 
could be due to the small amount of diffused Zn and a too 
short diffusion depth, as shown in Fig. 1. By contrast, the 
other two samples Nos. 2 and 3, treated at 600 and 700 °C, 
showed several dark modes in both TE and TM propagation, 
indicating optical waveguide formation was taking place. 

Table I shows a summary of the number of observed 
dark modes for samples Nos. 2 and 3, at the wavelengths 
provided by several laser lines. As can be seen, while sample 
No. 2 supports only TM propagation, and is monomode at 
\=488, 514, and 633 nm, sample No. 3, fabricated at 
700 °C, reveals several modes on all the wavelengths 
checked for both TE and TM propagation, indicating that 
both the ordinary and extraordinary refractive indices of the 
diffused layer have increased. The fact that the sample 
treated at 700 Q C reveals less number of modes for n e than 
for /i 0t in contrast with the observed modes in sample No. 2, 
is somewhat unexpected, taking into account that both 
samples have the same surface Zn concentration, and sample 
No, 3 is three times deeper. This behavior could be related to 
the difficulty of exciting all the propagating modes by prism- 
coupling, as Young et ai found 10 in waveguides fabricated 
by ZnO diffusion in congruent LiNb0 3 , and may be related 
to local variations in Li concentrations. 10 

Although the number of observed modes for TE propa- 
gation in sample Nos. 2 and 3 are too low to determine the 
exact index profile of the waveguides, this determination can 
be achieved for the ordinary index on samples treated at 
700 °C. Figure 2 presents these results, where we have plot- 
ted the increase in ordinary refractive index (as % to the 
unmodified substrate value) as a function of the penetration 
depth, obtained by fitting the dark mode spectra at \=488 
nm to a Gaussian profile A/i 0 — A/t 0|max exp(— x^iS). In 
order to show the goodness of this profile, we have also 
plotted in this figure (upper x axis) the calculated mode curve 
for this profile (dotted line) and the measured effective re- 
fractive indices (triangles). This fit gives a depth of d~2.7 
fjxn, in agreement with the diffusion depth extracted from 
EDAX data, and a value of An 0>max ^ 1.4%. Combining 
this value with the Zn concentration near the surface, a value 
of An 0 (%)/iZn]^030 is obtained. This ratio is about 
twice the value earlier obtained for ZnO:LiNb0 3 
waveguides, 10 the reason for this difference could be due to 
the different wavelength used for measuring the refractive 
indices (\=633 nm in referenced work, X=488 nm in this 
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FIG. 2. Gaussian refractive index profile at 4S8 nm for n<y (continuous line) 
of a Zn-diffused LiNbOj waveguide fabricated at 700 °C during 4 h. Dotted 
line represents the mode curve calculated for this profile, and triangles are 
the measured effective refractive indices (upper x axis). 

work) and the uncertainties for measuring the refractive in- 
dex prattle and the dopant concentration. 

The satisfactory results reported here are clearly at vari- 
ance with previous unsuccessful reports of waveguides pro- 
duced by vapor Zn diffusion in LiKb0 3 , a where the authors 
did not find evidence of waveguide formation by means of 
the dark mode characterization technique. To explain this 
failure one can possibly mention the use of too high a tem- 
perature for the diffusion, or the different method used to 
carry out the diffusion. 

It is worth noting that at the relatively low temperature at 
which waveguide formation was observed (600-700 °C) that 
any outdifrusion of Li ions is precluded. Thus, one of the 
drawbacks pointed out in Ti-diffused LiNb0 3 waveguides, 
where the diffusion is performed at 950-1200 °C, is over- 
come. Also, the low temperature diffusion is desirable to 
fabricate waveguide lasers in doped LiNbOj, thus avoiding 
the possibility of outdifrusion of the dopants. 

It has been established that addition of ~6% of Zn to the 
melt dramatically reduces the photorefractive effect in bulk 
LiNb0 3 . 12 Therefore, the use of Zn as an element to carry 
out the diffusion could help to avoid the photorefractive 
damage on these waveguides, as it has already been shown in 
Zn waveguides. 13 This element is an even better choice than 
Mg, because it does not modify the spectroscopic properties 
of doubly doped LiNbO s :RE:Zn. 14 

In conclusion, we have shown an alternative method to 
fabricate optical waveguides in LiNb0 3 by diffusing Zn from 
the vapor phase. Using metallic Zn for diffusion has the ad- 
vantage of the low vapor pressure of the Zn and the possi- 
bility of eliminating photorefractive damage. Good results 
have been obtained with diffusion carried out at temperatures 



as low as 600-700 °C, and the diffusion coefficient for Zn at 
700 °C is comparable to that obtained from Ti diffusion at 
1 1 00 °C, This low temperature processing is significant in 
that it slows down outdiffusion of Li and other elements in 
doped LiNb0 3 , and avoids perturbation of ferroelectric do- 
main distributions. These advantages may counteract the use 
of an experimental apparatus somewhat more complex than 
the traditional ZnO diffusion. Let us finally indicate that it 
has been recently demonstrated the fabrication of channel 
waveguides in LiNb0 3 by Zn diffusion from vapor phase, 15 
and the same method could be used for LiNb0 3 . The method 
here reported offers a versatile way to tailor the parameters 
of the guide by independently controlling temperature, time, 
and vapor pressure, although a precise prediction of the re- 
lationship between processing and waveguide parameters for 
engineering applications needs further characterization. 

Future work includes measurements of the photorefrac- 
tive damage threshold and losses of the produced 
waveguides, and investigating the possibility of using the 
method here described to extend the fabrication of optical 
waveguides to other insulator materials. 
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Preservation of Periodically Poled Structures in Zn-Diffused LiNb0 3 Waveguides 
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In this work, the formation of optical waveguides in periodically poled lithium niobate (PPLN) substrates by Zn-diffusion 
from vapour phase is studied. It has been found that the diffusion of the Zn ions is homogeneous in the PPLN substrates 
and the optical waveguides can be fabricated with the same parameters used in single-domain lithium niobate. The periodic 
ferroelectric domain pattern initially present in the substrate is preserved in the waveguides. 
KEYWORDS: lithium niobate, PPLN, waveguide. Zn-diffusion. microdomain 



l f Introduction 

Periodically poled (PP) non-linear crystals offer the pos- 
sibility of obtaining compact coherent light sources via opti- 
cal frequency conversion by quasi-phase matching (QPM). l,2) 
This procedure has been demonstrated using ferroelectric 
crystals such as KTiOP0 4 (KTP), LiB 3 0 5 (LBO) ] > 3 > and 
particularly in LiNb0 3 (LN) thanks to the large effective 
non-linear coefficient of periodically poled lithium niobate 
(PPLN). 3 > 

Taking advantage of the high power density available in 
wave-guiding structures, several nonlinear waveguide devices 
have been recently reported,^ in which the periodic modu- 
lation of the effective non-linear coefficient is combined with 
guided-wave geometries. 

LiNb03 nonlinear waveguides incorporating periodic do- 
mains have been fabricated using different techniques,* -12 ^ 
but these studies have identified difficulties associated with 
domain switching caused by the waveguide fabrication pro- 
cess 4, 13) at high temperatures. 

An alternative method to fabricate LiNb0 3 waveguides 
at relatively low temperatures has been recently demon- 
strated. 14) In this technique, Zn ions from a vapour bath are 
exchanged and then diffused using temperatures (700-900°C) 
substantially lower than the Curie temperature and therefore 
offers the prospect of preserving periodically poled structures 
previously defined, with potential advantages in integrated 
optics technology 

In this work it is demonstrated: (i) the feasibility of fab- 
rication of waveguides in periodically poled lithium niobate 
substrates by Zn diffusion, and (ii) the initial periodic ferro- 
electric domain pattern present in the substrate is preserved in 
the waveguides. 

2. Experimental Results and Discussion 

Figure 1 shows an optical microscopy image of the lithium 
niobate substrate used in this work. The substrate has been 
cut (z-cut) from a crystal boule grown in our laboratory by the 
Czochralski method using automatic diameter control system 
by weighing crucible system. 15 ^ 

The periodically poled pattern observed in the figure cor- 
responds to an area of around 2 x 2 mm 2 with a modula- 
tion period of 35 jtxm approximately This periodic structure 
was induced during the crystal growth by using the appro- 
priate growth conditions. 16) The ferroelectic domain pattern 
was revealed by etching using a diluted solution of HF:HN0 3 
(1:2 by vol) at room temperature. 4) This is subsequently pol- 
ished to a flat surface with optical quality and then the planar 




Fig. L. Optical microscopy photograph showing the periodic ferroelectric 
domain pattern in a z-cut lithium niobate substrate. 

waveguide was fabricated by Zn diffusion from the vapour 
phase l7 > at a temperature of 800°C. 

After diffusion the sample surface appeared clean and free 
of damage, allowing the optical characterisation by the prism 
coupling technique. l7) This has been performed using a lin- 
early polarised He-Ne laser (X = 632.8 nm) and a rutile 
prism. The assembly is mounted in a rotatory stage linked to 
a PC which allows the automatic determination of the sample 
reflectivity, with an angular resolution of 0.001 °. This proce- 
dure gives a precision of 1 x 10" 5 in the determination of the 
effective refractive indices of the modes. 17 ' 1S ) 

Figure 2 (continuous line) shows the reflectivity spectra of 
the waveguide taken at 632.8 nm after the Zn-diffusion pro- 
cess. As it can be observed, Zn-diffusion increases both the 
ordinary and extraordinary indices allowing TM and TE prop- 
agation (Figs. 2(a) and 2(b), respectively). The reflectivity 
curves indicate that the waveguide has a monomode char- 
acter for both propagations, in accordance with previous re- 
sults in single-domain LN substrates. I7) These results indicate 
that waveguide fabrication in PPLN is similar that in mono- 
domain LN substrates. 

After the initial optical characterisation stage, the waveg- 
uide has been carefully, etched by immersion in a diluted 
HF:HN0 3 solution (1 :2 by vol) at room temperature, and in- 
spected regularly at short time intervals (15 riiin.) in order to 
stop the etching immediately that the domain structure was 
observed. It was found that after 4 hours the ferroelectric do- 
main pattern is revealed sufficiently to identify the positive 
and negative domains. 
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Fig. 2. Reflectivity spectra of the optical waveguide taken at 632.8 nm 
(HeNe laser) for (a) TM and (b) TE propagation. Continuous and dot- 
ted lines represent the spectra before and after the etching process, respec- 
tively. 




Fig, 3, SEM micrograph of the PPLN structure after waveguide formation 
by Zn-diffusion. 



The waveguide was observed in a Scanning Electron Mi- 
croscope (SEM) to asses the quality of the ferroelectric do- 
main structures. Figure 3 shows a SEM micrograph cor- 
responding to the Zn-diffused PPLN where it is possible 
to recognise the alternate ferroelectric domains, which have 
been preserved with no changes in comparison with the struc- 
ture before waveguide preparation {Fig. 1). 

Nevertheless, it would be necessary to confirm that this 
structure corresponds to the waveguide and that this has not 
been removed by the etching process (otherwise the observed 
PP structure could correspond to the underlying undiffused 
region). This fact is confirmed by two procedures: First, En- 
ergy Dispersive X-Ray Analysis (EDX) measurements have 
been performed after the etching process in order to confirm 
the presence of Zn ions. This analysis indicates a concentra- 
tion of Zn ions close to 5.0mol% related to Nb ions, which 
is similar to the concentration in the waveguide prior to etch- 



ing. This concentration is the same in the positive and in the 
negative domains, indicating that the metal diffusion is homo- 
geneous in the PPLN substrate. 

Further information is obtained again by optical character- 
isation of the sample after the etching process. In fact the 
final confirmation of the waveguide presence after the etch- 
ing process is obtained by coupling the light to the sample 
using again the prism coupling technique. Figure 2 (dotted 
line) shows the reflectivity spectra after the etching process. 
It can be observed that the waveguide remains practically sim- 
ilar than before the etching process, which confirms that the 
waveguide parameters (depth and index change) are the same 
than prior to the etching and therefore the observed periodic 
structure corresponds to the waveguide region. 

So, in conclusion, it has been demonstrated that it is possi- 
ble the fabrication of optical waveguides in PPLN substrates 
by Zn-diffusion from the vapour phase, using the same dif- 
fusion parameters than in mono-domain LN substrates. Zn 
ions are incorporated homegeneously to the positive and neg- 
ative domains of the PP structure. The initial PP ferroelectric 
doimain pattern has been preserved in the waveguide by the 
Zn-diffusion process. These results are promissing for the 
fabrication of integrated OPO's and other frequency conver- 
sion devices based on lithium ntobate. 
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In this work, second harmonic generation (SHG) in nonlinear channel waveguides by quasi phase 
matching is reported. The photorefractive damage resistant waveguides were fabricated by 
Zn-diffusion technique using periodic poled lithium niobate wafers (PPLN). The SH wave, at 
397.8 nm, was generated by using a CW diode laser with emission wavelength at 795.6 nm and 
1 W peak power was used as fundamental wave. The nonlinear waveguides exhibit a stable conver- 
sion even at room temperature with no sign of deterioration induced by the optical damage. 

Introduction Due to the large effective nonlinear coefficient of periodically poled 
lithium niobate (PPLN) this crystal provides an attractive route for the realisation of 
compact coherent light sources via optical frequency conversion by quasi-phase match- 
ing (QPM) [1, 2], Combining the periodic modulation of the effective nonlinear coeffi- 
cient with the high power density in wave-guiding structures, several nonlinear wave- 
guide devices have been reported [4-7]. 

The method we have used to fabricate LiNb0 3 waveguides is a two-step Zn-diffusion 
[8], This technique permits to use PPLN substrates in which the periodic modulation is 
induced during the crystal-growth process. It is also important to consider that Zn-diffu- 
sion from the vapour phase at low temperatures produces low-loss optical waveguides 
with high photorefractive damage resistance [9-10]. 

In this work SHG, at room temperature, in the blue spectral region has been detected 
for fundamental wave in the near-infrared range using a CW diode laser to investigate 
the optical frequency conversion by QFM in Zn-diffused LiNbOs waveguides. 

Experimental Procedure Congruent lithium niobate crystals ([Li]/[Nb] = 0-945) peri- 
odically structured were grown along the jc-axis by the Czochralski method with auto- 
matic diameter control by weighting the crucible. This technique allows to induce the 
formation of periodic domain structures during the crystal growth using the appropriate 
growth conditions as reported elsewhere [11], 

The crystal boule was X-ray orientated and then wafers with the selected crystal 
orientation were cut. These wafers were polished to a flat surface with optical quality 
prior to perform the channel waveguide fabrication. 

Ultraviolet photolitographic techniques, adapted from microelectronic technology, 
have been used to define the channels in which the metal diffusion takes place [12]. 
Once the channels are defined, the waveguides were fabricated by a two-step Zn diffu- 
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Fig. 1, Experimental setup to obtain the 
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filter sion technique [8], In the first step 
the sample was placed in a Zn 
atmosphere, with pressure controlled 
by an Ar buffer, and heated during 
2h at 550 °C As a consequence of 

* electrometer this process a Zn-rich layer is 

formed in the surface of the sample. 
In the second step, this layer is diffused by annealing the sample in open air during 4 h 
at 850 °C This method produces low-loss optical waveguides preserving the periodic 
structures induced during the crystal growth [12]. 

In order to obtain the second harmonic signal the following experimental setup 
(Hg. 1) was arranged. A pig-tailed AlGaAs CW diode laser, with emission peaking at 
795.6 nm t was coupled to the channel waveguides by using a 40 x microscope objective. 
The output light from the waveguide was collected using a 20 x microscope objective 
and directed to the entrance slit of a monochromator (ARC SpectraPro 500-1) where 
both signals, excitation and SHG, were discriminated and recorded by online computer- 
controlled detection system using an EM1-9558QB photomultiplier tube. To avoid con- 
tamination of the SHG signal from the intense laser excitation a cut filter (LS-450-S- 
1075-CORION) was placed just in front of the monochromator slit. 

Results and Discussion Figure 2 shows an image of the waveguides by reflective light 
optical microscopy. Two different waveguides, 10 and 20 u.m wide, can be observed. The 
microscope has been focused onto the crystal surface in order to show the periodic 
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Hg. 2. Image of a PPLN crystal observed by an optical microscope with reflective light 

Fig. 3. Blue emission band corresponding lo the SH wave generated in the Zn-diffuscd nonlinear 
waveguides for a fundamental wave centered at 795.6 nm 
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Fig. 4 (online colour). Quadratic 
dependence of the SH intensity ver- 
sus the power of the fundamental 
beam 



0.01 



Fundamental power (mW) 



domains perpendicular to the channel waveguides, whose walls appear then diffused/ 
unfocused. The periodicity of the domains is 5 urn, which corresponds to a high effi- 
ciency of the SHG at a wavelength X = 397.8 nm [13]. 

In Fig. 3 the second harmonic signal can be observed. The emission of the blue light 
has been associated to TM— >TE conversion by QPM, in accordance with the phase- 
matching condition [13] taking into account the periodicity of the substrates. The funda- 
mental wave was polarised along the ordinary direction, in order to produce an interac- 
tion with the di2 coefficient, giving a strong polarised character along the extraordinary 
direction for the second harmonic wave. 

Figure 4 shows the dependence of the blue light intensity on that of the fundamental 
beam. The least square Gtting shows that it follows a quadratic dependence, as ex- 
pected for the nonlinear conversion mechanism. 

For a coupled pump power of 1,28 rnW, which corresponds to a value slightly lower 
than the maximum coupled power attained within the waveguide, the power obtained 
in the blue was 78 nW representing a conversion efficiency of ~10~ 4 . This value is 
similar to that previously reported for Ti-diffused LiNbOj waveguides [3], and lower 
than theoretical estimations, which could be attributed to a non-uniform periodicity of 
the domain that implies that the effective interaction length is reduced [13]. 

At variance with Ti-diffused waveguides [3] the infrared to blue light conversion ob- 
tained in this work using Zn-diffused waveguides, was operated at room temperature 
with a stable output and no sign of deterioration in the SHG. Also, as it has been 
already mentioned, the power dependence was quadratic during the whole range of 
excitation power used in these experiments, and no sign of photorefractive damage. 
These results confirm the potential advantages of the fabrication of channel waveguides 
by Zn diffusion, although the efficiency should be increased by improving homogeneity 
and a matching period of the domains. 

Conclusions It has been demonstrated that the periodic structures induced during the 
crystal growth in LiNbC>3 allow SHG under AlGaAs pumping at A = 795.6 nm, using 
Zn-diffused channel waveguides. These experimental results prove that the nonlinear 
characteristics are preserved by Zn diffusion and simultaneously the photorefractive 
damage is reduced. 
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ABSTRACT In this work second-harmonic generation by quasi- 
phase matching (QPM) in Zn-diffused periodically poled 
lithium niobate channel waveguides is presented. A stable 
TM-*TE conversion by QPM has been Found The results are 
in good accordance with theoretical estimations obtained by the 
phase-matching condition, either for the polarisation character 
of the second-harmonic wave as well as for the spectral range, 
taking into account the periodicity of the domains. 

PACS 42,82.-m; 42.65.Ky; 77.S4.-s 



1 Introduction 

Optical frequency conversion by quasi-phase 
matching (QPM) has been demonstrated using different fer- 
roelectric crystals [1-3]. Particularly interesting is lithium 
niobate (LiNbO*) due to the large effective non-linear co- 
efficients of periodically poled lithium niobate (PPLN) [1]. 
Several non-linear devices based on this material, in which 
the periodic modulation is combined with the high power 
densities available in waveguiding structures* have been 
reported [4-7]. 

The potential advantages of Zn diffusion have been pre- 
viously reported [8,9] as a fabrication technique to obtain 
LiNbO* waveguides. It has also been recently reported that Zn 
diffusion from the vapour phase at low temperatures produces 
low-loss waveguides in this material, and preserves the ini- 
tial domain pattern of the substrate [10, 1 1]. Therefore, with 
this technique it is possible to use PPLN substrates in which 
the periodic modulation is induced during the crystal-growth 
process. 

In this work the non-linear character of Zn-diffused pe- 
riodically poled LiNb03 channel waveguides is explored. 
A CW Ti:sapphire laser has been used to investigate the op- 
tical frequency conversion by QPM. Second-harmonic gener- 
ation (SHG), at room temperature, in the blue spectral region 
has been detected for fundamental waves in the near-infrared 
range, A good accordance between experimental data and the- 
oretical estimations has been found. 
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2 Experimental procedure 

Periodically poled lithium niobate crystals have 
been grown along the x-axis in an air atmosphere with 
automatic diameter control by a crucible weighing sys- 
tem [12], The initial melts, containing the congruent LiNbCb 
([Li]/[Nb] = 0.945) T have been doped with erbium and ytter- 
bium oxides with a purity grade of 99.99%. The concentration 
of both dopants in the initial melts was 0.5 mol%. This tech- 
nique allows us to induce the formation of periodic domain 
structures during the crystal growth by using the appropriate 
growth conditions [13- 15]. 

After cutting and polishing to optical grade, the ferroelec- 
tric domain pattern of a K-cut substrate from one of the PPLN 
crystals has been revealed by chemical etching in a diluted so- 
lution of HFiHNO} (1:2 by volume) at room temperature for 
1 h. Scanning electron microscope (SEM) measurements have 
been done to characterise the ferroelectric domains. 

The wafer was again polished lo a flat surface with optical 
quality in order to perform the channel-waveguide fabrica- 
tion. Ultraviolet photolithographic techniques, adapted from 
microelectronic technology, have been used to define the 
channels in which the metal diffusion takes place [16], The 
Zn diffusion was performed following the two-step proced- 
ure reported in [10]. In the first step the sample was placed in 
a Zn atmosphere at 550 °C for 2 h, and a second step of diffu- 
sion at 850 °Cfor4 h was selected in order to produce low-loss 
waveguides. 

A CW Tiisapphire laser, with a tunabiliry range from 
720 nm to 870 nrn, was used to perform the characterisation of 
the SHG in the Zn-diffused PPLN channel waveguides. The 
intensity of the second-harmonic wave was detected by using 
an EMI-9558QB photomultiplier tube. 

3 Results and discussion 

The ferroelectric domain pattern of the substrate 
was revealed by chemical etching and then analysed by using 
a scanning electron microscope. The SEM measurements 
confirm that different periodically poled (PP) regions, about 
1 mm in length and of different periods, were generated dur- 
ing the crystal-growth process. Figure la and b show the 
micrographs corresponding to two regions with periods of 
5 |xm and 10 pm respectively. Although there are different 
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FIGURE 1 SEM micrographs of PP regions with a S-[un and b JO-iim 
periods 



regions with periods ranging from 5 ixrn to 16 nm> most of 
these periodically poled regions have a period value around 
6—8 jxm. 

The initial stage in the fabrication of the channel wave- 
guides is the definition of the mask that stops the Zn diffusion. 
The motifs of aCr 2 03 commercial mask have been transferred 
to a Si02 film, previously deposited over the substrate, by UV 
photolithographic techniques [16]. 

The perpendicularity between the channels in the SK>2 
stopping mask and the PP domains in the LiNbCb substrate 
has been checked, prior to the Zn diffusion, by diffraction 
techniques. A helium-neon laser (632.8 nm) was used to illu- 
minate the sample surface that contains the channels; Fig, 2a 
shows one of the diffraction patterns obtained. The diffrac- 
tion in the vertical line corresponds to the periodic series 
of channels in the Si0 2 film, while the pattern in the ho- 
rizontal line is induced by the presence of a periodic modu- 
lation in the ferroelectric domains, showing that a good per- 
pendicularity between channels and PP domains has been 
achieved. 

In Fig. 2b a SEM micrograph of the Zn-diffused channels 
(with dimensions 1-cm long and 1-, 2-, 4-, 6-, 10-and 20-u.m 
wide), before removing the SiCh stopping mask, is presented. 
The channel waveguides are parallel to the *-axis (the c-axis 
perpendicular to the channels). 




FIGURE 2 a Diffraction pattern at 632.& run showing the perpendicularity 
achieved between the channels in the SiOj stopping mask and the PP do 
mains: vertical and horizontal patterns respectively, b SEM micrograph of 
Zn-diffused waveguides 2-, 4 % 6-, 10-and 20- tun widths) 



The non-linear character of the Zn-diffused channel wave- 
guides has been explored in continuous-wave mode at room 
temperature by using a CW Ti:sapphire laser. The funda- 
mental wave (o>), polarised along the ordinary direction (in 
order to produce an interaction with the dyi coefficient), was 
coupled into the waveguides by using a microscope objec- 
tive. A second objective, at the other end, was used to collect 
the second-harmonic wave (2co). A strong polarised charac- 
ter (higher than 90%) has been found along the extraordinary 
direction for the 2o> wave, indicating that the QPM involves 
a TM— >-TE conversion. 

Figure 3 shows -the phase-matching curve for the TMh- 
TE conversion measured in the 2-uju channel waveguide with 
1-cm length (containing non-uniform PP regions, as has been 
mentioned previously). As can be seen, the phase-matching 
condition is satisfied for fundamental wavelengths higher than 
770 nm, the 2a> wave having its maximum intensity for ex- 
citation at around 822 nm The inset in this figure shows the 
expected quadratic dependence of the SH intensity on the in- 
tensity of the fundamental wave. 

Figure 4 shows the calculated relationship between the fer- 
roelectric domain period {A) and the fundamental wavelength 
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wavelengths close to 820 nm. From (1), it can be deduced 
that the periods needed for this conversion should be close 
to 7 pjTi. This value is also in good accordance with the ini- 
tial SEM analysts of the substrate (mentioned above), which 
reveals the dominant period value (around 6-8 \im) in the pe- 
riodically poled regions. 

Preliminary measurements related to the photorefractive 
damage indicate that, at the incident power level used in 
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FIGURE 3 Phase-matching curve at room temperature measured in the 
2-pliti channel waveguide. The inset shows the dependence of the SH wave 
on the intensity of the fundamental wave 
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FIGURE 4 Calculated relationship, in the case of QPM by TM-*TE con- 
version, between the ferroelectric domain period and the fundamental wave- 
length 



(Xa,) at first order for QPM by TM- 
given by [3]: 



A ^ 



2(/i t (2co)-n 0 (o>)] 



*~TE conversion, which is 



(1) 



these experiments (around 0. 1 W), the fundamental mode pat* 
tern and the SH intensity remain stable, although a quantita- 
tive estimation of the optical damage threshold has not been 
obtained. Detailed evaluation of this figure is currently in 
progress, 

4 Conclusions 

The feasibility of Zn diffusion from the vapour 
phase to produce non- linear waveguides in lithium mo 
bate has been demonstrated. Optical frequency conversion 
by quasi -phase matching has been obtained in Zn-diffused 
LiNb03 channel waveguides via TM-*TE conversion. The 
SH wave exhibits a broad phase-matching curve in conse- 
quence of variations in the period associated with the PP 
ferroelectric domains. 

Further investigations in order to precisely control the do- 
main period during crystal growth are needed in order to ob- 
tain a higher efficiency in blue-light generation. 
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Abstract 

Single crystal samples of lithium niobate have been spatially patterned with photoresist, and subsequently domain 
inverted using electric field poling, to produce a range of two dimensional spatial domain structures. Differential etching has 
subsequently been carried oul using mixtures of hydrofluoric and nitric acids, at a range of temperatures between room 
temperature and the boiling point. The structures produced show very smooth, well defined, deep features, which have a 
range of applications in optical ridge waveguides, alignment structures, V-grooves, and micro-tips. Details are given of the 
fabrication procedures, and examples of structures are shown. © 1998 Elsevier Science B.V. All rights reserved, 

Keywords: Lilhium niobate: Differential etching; Ferroelectric domains 



1. Introduction 

Lithium niobale is a material of considerable in- 
terest to the optical, laser and communications indus- 
try, due to the lar^e values of nonlinear optical, 
electro-optic, piezoelectric and acousio-optical coef- 
ficients and figures of merit that are available [1], 
Current uses span the laser area, where fabrication of 
periodically poled material has achieved record effi- 
ciencies in quasi-phased matched nonlinear interac- 
tions [2], to the microwave communications sector 
for use in surface acoustic wave delay lines and 
filters [3]. In the telecommunications field, it is 
widely used for integrated optical modulators where 
switching rates of 10 GHz and above are currently 
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required for wavelength division multiplexed fibre 
optical networks. Finally, although regarded as a 
problem in many areas, it continues to maintain 
interest in the photorefractive community for optical 
storage, holographic and phase conjugate applica- 
tions. 

A number of techniques have been reported for 
machining and microstructuring in lithium niobate. 
These include reactive ion etching (RIE) [4], reactive 
ion beam etching (R1BE) [5], laser ablation, both 
above and below bandgap [6], and laser assisted 
chemical etching using chlorine [7] and KF [8], Wet 
etching has also been used for many years to reveal 
domain and defect structures in lithium niobate [91 
but not so far to our knowledge for intentionally 
producing ordered rnicrostructures subsequent to 
lithographically defined, spatially patterned, domain 
inversion. Such niicrostrucUiring has obvious appli- 
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cations for the fields of optics and optoelectronics: 
the ability to spatially define very smooth and pre- 
cise structures for lateral guidance (waveguides), and 
multiple arrays (photonic bandgap devices) will have 
immediate impact. 

The most common wet etchant reported has been 
a mixture of hydrofluoric (HF) and nitric acids 
(HN0 3 X which attacks the negative z face (-z) at a 
rate that is appreciably higher than the positive z 
face ( + z). Additionally, defect sites in the crystal 
lattice may also experience higher etching rates. This 
last property has been utilised recently to produce 
ridge waveguides for improved electrical-to-optical 
coupling in integrated optical interferometer devices 
[10]. Many such techniques require the use of a prior 
lithographic spatial defining step, which has implica- 
tions in the ultimate quality of structures obtainable 
due to mask fabrication limitations. In our technique, 
although domain inversion still requires a litho- 
graphic procedure, involving photoresist inhibition of 
the electric field required for domain inversion, the 
latter process is inherently a manipulation of the 
crystal local structure at the atomic, or unit cell level. 
Subsequent differentia] etching can therefore pro- 
duce structures that reflect this level of perfection, 
and extremely smooth features are therefore generic 
to this technique. 

We report here the combination of spatially selec- 
tive poling of lithium niobate, and subsequent differ- 
ential etching, to produce structures up to 150 jLtm 
deep with very smooth walls, as imaged using scan- 
ning electron microscopy. The results obtained sug- 
gest that this technique may be very suitable for 
production of high quality ridge waveguides, V- 
grooves and other alignment aids, electro-optically 
active tips for atomic force microscopy applications, 
and arrays of highly uniform holes, for filtration and 
sensing purposes, and columns for photonic band 
gap structures. 



2, Materials and methods 

The samples of lithium niobate had dimensions of 
15 mm X 15 mm and were either 300 or 500 /im 
thick. All material was supplied by Koto Electric 
(Japan) in the form of z cut, 76 mm diameter 
wafers, which were cut to size in house. The — z 



face was marked to show the y direction, (Nomenc- 
lature of the coordinate system is as described in 
Ref. [1 1]). AH samples were subjected to a sequential 
multiple solvent cleaning process using ultrasonic 
agitation at 50°C. The samples were then dried with 
compressed air. Photoresist was spin coated onto the 
— z face, producing a LI /xm thick film. The sam- 
ples were finally baked at 90°C for 30 min. 

The photoresist was patterned using a Karl Suss 
MA4 mask aligner and chrome on quartz masks. The 
patterns used to measure the temperature dependant 
etch rate consisted of gratings with periods between 
3 and 27 /Am, channels of width between 3 and 14 
jam, and arrays of square, triangular and round 
points with feature sizes of 5, 10 and 20 Aim. The 
channels and gratings were aligned with their grating 
vector parallel with the crystal! ographic x direction, 
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Fig. 1 . (a) Schematic of the poling apparatus used, (b) Schematic 
of domain structure within the LiNb0 3 sample before and after 
poling. 
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and features parallel with the y direction. Other test 
patterns used were 200 ju,rn diameter circles, ran- 
domly orientated Hnes of 100 fxm width, and grat- 
ings with periods between 5 and 18 Atm. These latter 



structures were etched only at room temperature. 
Once exposed the resist was developed for 50 s, 
washed in de-ionized water for 60 s and dried using 
dry compressed air. 
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The samples were poled by applying an external 
electric field of 22.5 kV/mm. Fig. la is a schematic 
of the apparatus used. Poling current was monitored, 
and controlled by varying the applied field. The 
expected displacement current, and hence total 
charge, was calculated in advance, and the poling 
process was terminated once the required degree of 
poling had been achieved. Liquid electrodes were 
used to ensure good electrical contact with the sam- 
ple. After poling all samples were again cleaned 
using acetone and water 

The etchant consisted of HF and HN0 3 acids in a 
1:2 ratio, where the HF was a 48% solution and the 
nitric acid was a 78% solution. The etchant was 
slirred constantly using a magnetic stirrer with an 
integral hotplaic. The temperature of the hotplate 
was kept constant and the etchant allowed to reach 
equilibrium temperature before initiating the etch 
process. The etching lime varied between 15 and 17 
h lor all samples. The etch depth was measured 
using a Tencor Alphastep profileorneter. Addition- 
ally a JEOL JSM 6400 Scanning Electron Micro- 
scope (SEM) was used for visual inspection of the 
siruciures produced. Some samples were also set in 




Fig. 3. Simultaneous view of both faces ( + z and — z\ and different etch results. The edge of the sample appears rough, as this has been 
snapped prior to SEM viewing. 



low viscosity resin and then cut and polished in a 
plane normal to the z face to allow measurement of 
the etch depth from the SEM pictures obtained. 



3. Results and discussion 

Before poling, the lithium niobate wafers were all 
oriented single domain samples. After poling, sam- 
ples had areas of induced domain reversal, as is 
illustrated schematically in Fig. lb. Examination of 
the resultant etched structures revealed that opposing 
crystal z faces displayed quite different topology. 
The features on the photoresist patterned —z face 
appear to faithfully replicate the areas of induced 
domain inversion. The features on the opposing z 
face (the +z face), however, show some evidence of 
domain spreading, as revealed by the subsequent 
etching procedure. The underlying crystal 3 m sym- 
metry leads to features with characteristic triangular, 
and hexagonal shapes, and, as seen in the example of 
Fig. 2, an initially circular photoresist pattern will 
produce an equivalent circular etched structure on 



>: <XP 4256024A I > 



250 



/.£. Barry et al. / Materials Letters 37 (1998) 246-254 



the — z face, but a hexagonal structure on the etched 
+z face. 

Fig. 3 also shows such complementary structures, 
and illustrates the extent to which the etch process 
can be utilised to construct alignment aids, such as 
* V'-grooves. It should be noted that the exact profile 
is not that of a * V* in this case, but rather a tapered 
'V\ The reasons for this exact structure are currently 
under further investigation, and are being studied as 
a function of both etch ant concentrations and tem- 
perature. Under some conditions, very small (~ 100 
nm) holes appear around features revealed on the — z 
face. These are illustrated in Fig. 4 t and are also 
under active study for possible uses in filtering, 
sieving, etc., but with the additional benefit of being 
in an active optical material. We are currently con- 
sidering the modifications that can be expected in the 
nonlinear and electro-optic response of lithium mo- 
bate, by the inclusion of nanoscale particles of met- 
als, or dielectrics. As a general observation, all struc- 
tures produced via etching, and revealed on the — z 
face display smooth walls, and vertical features. 
Under SEM examination, the smoothness persists at 
the highest resolution available with our SEM (~5 



nm). The implication here is that crystal planes are 
being revealed via this technique and we are cur- 
rently measuring individual facet orientations. Fur- 
ther studies will involve X-ray diffraction and topog- 
raphy analysis. 

As etching proceeds, the facets produced appear 
to result in self- termination, or a much reduced rate 
of etching. This is of most use when alignment 
grooves are required, and an example of this is 
shown in Fig. 5. The facetting does not occur at the 
same depth in all cases however, and in Fig. 6 we 
show an example of a deep vertical walled ridge, 
suitable for waveguide applications, in which the 
etch has proceeded to a depth of ~ 20 /im, without 
any facets appearing. This behaviour is critical to the 
use of this technique for fabricating high aspect ratio 
structures. Finally we note the appearance of most of 
the etched samples so far in the near surface region. 
It is known that the domain inversion process is 
influenced by defects, and such structures as sub- 
surface polishing damage. The first few micrometers 
of all etched samples shows a 'banded* appearance, 
the detailed structure of which may vary slightly, but 
which appears to be common to all etched samples. 




Fig. 4, Holes produced around edges of etched domain walls. 
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Fig. 5. Etched structure showing smooth walls, and V-groove features. 




Fig. 6. Example of ridge structure suitable for waveguide applications. Note the smooth side- walls produced and clear evidence of surface 
'banded' appearance. 
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Fig. 6 shows this particularly clearly. The surface 
smoothness seems unaffected however, and we are 
further investigating the origin of the contrast ob- 
served. A final example of the precision that can be 
achieved by this method is shown in Fig. 7. An 
isolated tip has been produced which may well prove 
useful in scanning microscopy applications, such as 
near field imaging. 

Fig. 8 shows the results of etch rate obtained as a 
function of etch temperature. The mechanism by 
which etching proceeds faster at elevated tempera- 
tures can be considered from two perspectives. 
Firstly, such etching processes are expected to follow 



the normal Arrhenius law, (rate proportional to 
>lexp[ — E a /RT]\ which would yield an exponential 
increase in etch rate. Secondly, however, there is 
another factor which might influence temperature 
dependent etch rates. Poling acts to displace both 
Li + and Nb + ions relative to the oxygen ions in the 
lattice. As this displacement is the key indicator of 
domain polarity, the exact positions of these ions 
must play an extremely important role in subsequent 
differential etch rates. Increased temperature will 
alter the metal ion positions within the lattice, and 
hence may well influence etch rate. Fig. 8 shows an 




2,806-3 



3.00e-3 
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Fig. 8* Arrhenius plot showing natural tog of eich-rale (measured 
in /im per hour), vs. 1 /absolute temperature. The straight line fit 
shows no departure from the normal exp( £" a / RT) behaviour. 

Arrhenius plot of In (etch -rate), measured in /xm per 
hour, vs. 1 /absolute temperature, from which the 
linear gradient indicates no obvious departure from 




Fig. 7. Single crystal tip, for possible scanning microscopy applications. 
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Fig. 9. Deep etching of trenches. The — 150 ftm depth demonstrates independence of facet angle on number of separate etching steps, and 
feasibility of deep structures. 



the expected Arrhenius law. Five temperatures from 
21 to 95°C were used, but at the highest value, the 
etch rate was sufficiently fast to completely etch 
through the sample, preventing a quantitative deter- 
mination of etch rate. Fig. 8 also indicates that at the 
highest temperature that can be used for boiling 
HF/HNO3, ft-hich is 110°C, an etch rate of -55 
ju.m per hour is predicted. Repeated immersion of the 
sample did not affeel the observed facet angles. The 
sample seen in Fig. 9 was immersed twice, and no 
change in etch behaviour was observed. 

In summary therefore, a process consisting of 
domain patterning by electric field poling followed 
by etching in a 1:2 mixture of HF and HN0 3 acids 
has revealed a range of micro- and nanoscale struc- 
tures. The surface quality appears extremely smooth 
and straight in directions parallel to the mirror planes 
of the crystal 3 m symmetry, with surface roughness 
less than the 5 nm resolution available to us. The 
technique has produced ridge waveguides, alignment 
grooves, sharp tips and regular arrays of similarly 
shaped and equally sized holes. We are currently 
evaluating the technique for producing low loss ridge 
waveguides, and have defined a lower confinement 



layer using ion-beam implantation. The tip structures 
are under consideration for electro-optic sampling, 
near-field imaging microscopy, and other associated 
techniques. 
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Fabrication of LiNb0 3 TE/TM Waveguides for 1.5 /tm Wavelength Band 
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Zn-d iff used waveguides can support both ordinary and extraordinary polarization modes and are expected to have high resis- 
tance to photorefractive damages. In this paper, fabrication of ZniLiNbOj channel waveguides for 1 .5 Mm wavelength band by 
diffusing metallic Zn in low-pressure atmosphere is reported. The waveguides were fabricated under various conditions. It was 

Single-mode TE/TM waveguides with — I clB/cm propagation were obtained. The process is simple and has good reproducibil- 
ity/ IDOi: KUI43/JJAP.4L4825) 

KEYWORDS: lithium niobate, waveguide, Zn-diffusion, integrated optics, nonlinear optics 



1. Introduction 

LiNbQj waveguide quasi-phase-matched (QPM) nonlinear- 
optic (NLO) wavelength conversion devices have been stud- 
ied for applications in optical communication systems, optical 
signal processing, etc. Annealed-proton-exchanged waveg- 
uides are widely used for such devices ,,2) because of their 
high resistance to photorefractive damage. However, they 
support only extraordinary guided modes, and they may suf- 
fer limitations in implementation of polarization-independent 
devices. Ti-dif fused waveguides, which can support both or- 
dinary and extraordinary guided modes, are also used. ? 4> But 
they involve notable photorefractive damage problem. The 
formation of unwanted domain inverted thin layer during the 
high temperature process for Ti diffusion may cause problems 
in fabrication of domain inverted gratings for QPM. 4) 

Zn diffusion in LiNbOj takes place at relatively low tem- 
perature, and it can provide Zn-diffused waveguides which 
can support both polarization modes and are expected to have 
high resistance to photorefractive damages. 5 * They are poten- 
tially suitable for applications to QPM -NLO devices/'* Wc re- 
ported recently Zn-diffuscd LiNbO;* waveguides for 1.5 //m 
wavelength band by diffusing from ZnO 71 and their applica- 
tion to QPM wavelength conversion devices. Si There is also a 
report on the fabrication of the waveguides in y-cut LiNbO* 
for 1.32 /xm wavelength by diffusing metallic Zn for appli- 
cations to electrooptic (EO) devices. 91 However, fabrication 
techniques of Zn-diffuscd waveguides have not been fully es- 
tablished yet* especially for z-cut LiNbO.v 

tn this study, we fabricated Zn-diffused LiNbO_i waveg- 
uides for 1.5 /urn wavelength band in z-cut LiNbOj by dif- 
fusing metallic Zn for applications to NLO devices. Pressure 
control during diffusion was crucial to obtain smooth sur- 
face after diffusion. We fabricated the waveguides under var- 
ious conditions and clarified the condition to obtain TE/TM 
waveguides with good quality for 1.5 fim wavelength band. 

2* Fabrication 

Diffusion of metallic Zn into congruent optical-grade 
LiNbOj was examined. Although we tried to deposit Zn di- 
rectly on the -f-z surface of LiNbO* crystal by thermal evapo- 
ration, the surface of the deposited Zn film was granular. In- 




(a) 



< > 



10 |xm 



(b) 



Fig. I. Optical microscope photographs of the surfaces after diffusion ;il 
89Cr*C for 30m»n in (a) atmospheric pressure and fb) l.OTorr. 



sertion of a thin Ni buffer layer was reported to allow deposi- 
tion of the Zn film with uniform surface. 9 * We first deposited 
Ni film of 5nm thickness, and then deposited Zn on it. The 
Zn/Ni film was diffused using an ordinary quartz-lube elec- 
tric furnace. The thermal diffusion in air of atmospheric pres- 
sure resulted in rough surface, as shown in Fig. 1(a). !l was 
probably caused by Zn-UNbOji reaction in the air. We tried 
the thermal diffusion in air of low-pressure atmosphere. 7 * The 
quartz lube was evacuated by a rotary pump. The pressure was 
monitored with a Pirani gauge and was controlled by a throt- 
tle valve, Wc found that the smooth surface, as shown in Fig. 
1(b), was obtained by diffusion in a low-pressure atmosphere 
of about 1 .0 Tom Thermal treatment in low-pressure acceler- 
ated Zn diffusion into LiNb0 3 . Similar tendency was reported 
for Zn-diffuscd waveguide fabrication in vapor phase. 10 * The 
pressure higher than KOTorr resulted in rougher surface, and 
the lower pressure caused Li 3 0 outdiffusion, 

Zn/Ni stripes of 4-8 fj.m width were patterned by lift-off 
technique, and diffused under various conditions for channel 
waveguide fabrication. The total thickness of Zn/Ni film, the 
diffusion temperature, the diffusion time ranged 1 00-200 nm T 
830-95O°C\ 30-120min, respectively. The pressure ranged 
from 1 .0 to 30Torr. After the thermal diffusion, the end facets 
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Table I. Result of the optical guiding test of the fabricated chunneJ waveguides. 



Fabrication conditions Waveguiding characteristics 
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Fig. 2. The surface of Zn-diffused channel waveguide. 



were polished for end-fire coupling. The waveguide length 
was approximately 10 mm. The typical surface morphology 
of the crystal after diffusion is shown in Fig. 2. This sample 
was fabricated by diffusion at 890°C for 30rnin in I.OTorr. 
Relatively smooth surface was obtained, though there was 
some residual thai might cause a little propagation loss. 

3, Characterization 

The waveguiding characteristics at 1.55 wavelength 
were examined. A beam from a hber-pigtailed laser diode 
was coupled into a waveguide using a 20 x lens. We eval- 
uated each waveguide by the full width at half-maximum 
(FWHM) mode size and throughput. The results are summa- 
rized in Tabic L In case of diffusion at 890*C\ the waveguide 
#3 supported both polarization modes with relatively high 
throughput. Fn cases of diffusion at 830°C (#1) and 860°C 
(#2), TE mode showed quite low throughput. The waveguide 
#4 fabricated by diffusion for 120 min had larger mode size; 
they didn't have so good confinement as the waveguide #3 
of 30 min diffusion. It seems to be because the diffusion was 
too deep. It was seen by comparing #3 and #5 that the Zn/Ni 
thickness should be about I60nm. The waveguides diffused 
from Zn/Ni stripe of 8 fim width were multi-mode waveg- 
uides for TM mode. To make single-mode waveguides for 
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Fig. 3. The dependence of throughput on diffusion temperature. 



both polarization modes, Zn/Ni stripe width should be about 
5 jj.m. In case that the width is 5 /xm, Zn/Ni stripe thickness 
needs to be about 200 nm. In addition, diffusion temperature 
has to be high. Under such conditions, the diffusion in I.OTorr 
progresses too rapidly and resulted in poor optical confine- 
ment. To make the diffusion rate slower, the pressure should 
be increased to approximately 20Torr. 

Figure 3 shows dependence of throughput on diffusion tem- 
perature. The higher the diffusion temperature was, the higher 
throughput and the smaller throughput differences between 
TM and TE modes we obtained. The reason for the depen- 
dence is under study. It is seen from Fig. 3 thai the diffusion at 
higher temperature is favorable for the fabrication of waveg- 
uides with high throughput for both polarization modes. The 
best condition to fabricate single-mode channel waveguides 
with good confinement and high throughput for both polariza- 
tions is diffusion of Zn/Ni stripe with 5 yum width and 200 nm 
thickness at 950°C for 30 min in 20Torr air. 

Figure 4 shows the mode patterns in the waveguide fabri- 
cated by diffusion under the best condition. The FWHM mode 
sizes were 7.GM r n(depth) x 9.7 //m(width) 1 and 7.8 ^xm x 
9.7 /xm for TM and TE modes, respectively. The propa- 
gation loss was measured as ~~ \ dB/cm for both modes 
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(a) 




by the waveguide Fabry- Per<5t interference method using a 
temperature-tuned distributed feedback laser. 

A, Conclusion 

We studied fabrication of Zn-diffused LiNbO^ waveguides 
by diffusing metallic Zn. and we found that diffusion in low- 
pressure atmosphere is crucial to fabricate the waveguides. 
We fabricated channel waveguides under various conditions 
and examined the wavcgujdjng characteristics at 1.55 jam 




wavelength. Single-mode TE/TM waveguides for 1.5 ^m 
wavelength band with \ dB/crn propagation loss were ob- 
tained. The process is simple and has good reproducibility. 
Future work includes examination of the resistance of the 
Zn -diffused waveguides to photorefractive damages. More re- 
cently we applied these Zn-diffused waveguides for QPM 
wavelength conversion devices, and obtained encouraging re- 
sults. The results will be reported elsewhere. 
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Fi£, 4. Near field patterns of guided modes at I.S5/im for (a) TM mode 
and (b) T\i inodc. 
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Abstract 



In this paper, photorefractive crystals of Ce,Fe:LiNb0 3 are systematically studied. The crystals have been grown by the 
Czochralski method. Samples with various doping concentrations and oxidation/reduction treatment have been fabricated. 
Their photorefractive properties were experimentally investigated by using two-beam coupling. The results show that, 
compared with Ce:LiNb0 3 and Fe:LiNb0 3 , the photorefractive effect of Ce,Fe:LiNb0 3 crystal was enhanced, and the 
photorefractive efficiency depends on the dopant concentration, oxidation/reduction treatment, and light wavelength. The 
doping mechanism is also discussed here. ©2000 Published by Elsevier Science B.V. All rights reserved. 

Keywords: Photorefractive properties; Two-beam coupling; Ce>Fe:LiNbO s crystal 



1, Introduction 

Doped lithium niobate (LiNb0 3 ) crystals exhibit 
excellent photorefractive properties, and this attrac- 
tive feature makes them important for applications in 
optical computing, image processing, phase conjuga- 
tion, and laser hardening. Doping with transition 
metal or rare earth ions into crystals can affect 
efficiency, sensitivity, speed and spectral response of 
photorefractive effect. A great amount of studies on 
LiNb0 3 included crystals containing various dopants, 
such as Fe [1], Cr [2], Nd [31 Pr [4], Eu [51 Mg [6], 
Zn [71 Ce [8], and so on, Fe: LiNb0 3 crystal has 
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been selected as one of holographic storage media 
for industrial use. It is due to the outstanding advan- 
tages as follows: (1) it can be easily grown from 
melt by Czochralski method; (2) its Curie tempera- 
ture is very high (around 1200°C), so it is convenient 
during processing and applications; (3) its photore- 
fractive effects are strong; and (4) there are some 
vacancies in the lattice, which makes it easy to 
introduce some dopants to tailor its photorefractive 
properties. 

The effect of different dopants on the photorefrac- 
tive effect of crystal is different Two kinds of 
dopants were doped in LiNb0 3 in order to enhance 
its photorefractive effect, for example, the response 
time of Mg,Fe:LiNb0 3 is one magnitude faster than 
that of Fe:LiNb0 3 [91 Double-doped Ce,Fe:LiNb0 3 
crystal was grown in our lab in 1988 t and has been 
proved as an excellent photorefractive material by 
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the experimental results of two-beam coupling, 
four- wave mixing, and image processing [10-12]. In 
this paper, the growth and photorefractive properties 
are systematically reported, 

2. Crystal growth and sample preparation 

The starling materials used for crystal growth are 
ultra pure Li 2 C0 5 , Nb 2 O s , Ce0 2 and Fe 2 0 3 . The 
congruent composition of Li:Nb = 0.946 was se- 
lected as melt and crystal compositions* Ce ions 
were added into the melts with the concentration of 
0.1 rnol.%, and Fe ions were with varying concentra- 
tions of 0.02 and 0.05 mol.%. The appropriately 
weighed materials were thoroughly mixed and cal- 
cined at 700°C and then sintered at 1150°C for 2 h, 
respectively. 

Single crystals of Ce:LiNbO s , Fe:LiNb0 3 and 
Ce,Fe:LiNb0 3 were pulled along c axis from the 
melt in a platinum crucible by Czochralski method. 
The sizes of crystals as grown are ~ 30-35 mm in 
diameter and 30-40 mm in length. Typical technique 
parameters for crystal growth are listed in Table 1, 

All as-grown crystals show good symmetry and 
optical quality. The crystals were placed in a furnace 
where the temperature gradient is almost equal to 
zero for annealing and poling. After annealing at 
1200°C for 8 h, the crystals were polarized with 
electric field at 5 mA/cm 2 for 30 min at the same 
temperature. Then the crystals were cut into samples 
with size of 12 X 12 X 2 mm 3 (a X b X c\ 

As well-known, modest treatment of oxidation 
and reduction can change the valence state of the 
doped ions and the other defects* distribution, there- 
fore affect the photorefractive effect of the crystal. In 
our experimental, in order to avoid introducing other 
impurities during oxidation and reduction, the pow- 
ders of Nb 2 Os Li 2 CO s were chosen as the 



Table 1 

Growth conditions for Ce:LiNb0 3 , Fe:LiNb0 3 and Ce T Fe:LiNb0 3 
crystals 



Growth Temperature 


Growth 


Growth 


Crystal 


temperature gradient 


atmosphere 


rate 


rotation 








rate 


1250°C 40°C/cm 


air 


2-3 mm/h 


15-25 rpm 



Table 2 

Samples of Ce:LiNb0 3 , Fe:LiNb0 3 and Ce,Fe:LiNb0 3 crystals 



No. 


Ce0 2 
(mol.%) 


Fe 2 0 3 
(moJ.%) 


Treatment 


Color 


1 


0 


0 


as grown 


pale yellow 


2 


0.10 


0 


reduction 


pale brown 


3 


0 


0.05 


reduction 


red 


4 


0.10 


0.02 


as grown 


rose red 


5 


0.10 


0.02 


reduction 


heavy red 


6 


0.10 


0.05 


as grown 


red 


7 


0.10 


0.05 


oxidation 


pink 


8 


0.10 


0.05 


reduction 


brown red 


9 


0.10 


0.05 


heavy reduction 


red black 



oxidizing and reducing agents, respectively. The 
sample sheets were buried in the powdered Nb 2 0 5 
(or Li 2 C0 3 X The temperature of the oxidation was 
1150°C, the treatment times were 24 h. The tempera- 
ture of reduction and heavy reduction were 550°C, 
the treatment times were 24 and 36 h respectively. 
The treated samples were lapped and optically pol- 
ished for characterization (see Table 2), 



3. Optical absorption properties 

The absorption spectra were obtained with non- 
polarized light and shown in Fig. 1. There is a broad 
absorption band located at 460-550 nm, which is 
induced by Fe 2+ ions and oxygen vacancies V£ [131 
It is obviously that with more Fe ions in the crystal 
the absorption becomes heavier. 

The influences of oxidation and reduction are also 
clearly observed. The intensity of the absorption 




° 400 500 600 700 &00 
A (nm) 

Fig. 1. Absorption spectra of CeO 2 (0.10 moL%),Fe 2 O 3 (0.50 
moI.%):LiNb0 3 crystals: (a) No. 9 sample with heavy reduction; 
(b) No. 8 sample with reduction; (c) No. 6 sample as grown; and 
(d) No. 7 sample with oxidation. 
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band at 460-550 nm increases as reduction, and 
decreases as oxidation. That is resulted in the reac- 
tions that Fe 2+ ions and oxygen vacancies V£ are 
induced by reduction treatment, but eliminated by 
oxidation treatment. The effects of these treatments 
on photorefractive properties will be studied and 
discussed later. 




Fig. 2. Experimental set-up of two-beam coupling, BS: beam 
splitter, M: reflect mirror. 



4. Two-beam coupling: theory 

When two coherent laser beams incident into a 
p hot orefr active crystal y in wh i ch t he C axis was 
oriented to be in the incident plane and perpendicular 
to the bisector of the two beams, beam coupling 
occurs, and causes one beam (signal beam) to gain 
intensity at the expense of the other beam (pump 
beam). The exponential gain coefficient of two-beam 
coupling r is experimentally measured through the 
relation [14,15], 



1 K 1 * 



where L is the interaction length of the two beams in 
crystal (taken as the thickness of the sample), (/ s ) 
is the transmitted signal beam intensity with (without) 
coupling, and / p (/ p ) is the transmitted pump beam 
intensity with (without) coupling. When the intensity 
of pump beam is much larger than that of signal 
beam, / p :s> / s+ the two-beam coupling gain coeffi- 
cient r becomes independent of the pump beam 
intensity (i.e. fp~I p \ 



i /; 



The two-beam coupling gain coefficient T is related 
to the photorefractive index grating amplitude Sn 
(defined as one-half the peak-to-peak value) by [16] 



f—A-ir 5/i sin <J>/m A cos 0j 



(3) 



where 6 1 is the half angle between the beams inside 
the crystal, </> is the phase shift between the optical 
interference pattern and the photorefractive index 
grating, and m is the modulation depth of th e inci- 
dent optical interference pattern, m = 2^// p / s /(/ p 4- 



/ s ). The photorefractive index modulation 5/i is given 

by 

(4) 



Sn = n 3 r eff £ sc /2 



where E sc is the space-charge electric field, n is the 
refractive index, and r eff is the effective electro-optic 
coefficient, The photorefractive space-charge field 
E 5C is given by [17] 

E sc =m(k B T/e)[K/(l 4- ( K/K 0 ) 2 ] f(^)cos2^ 

(5) 

where 20, is the internal full crossing angle of the 
optical beams and k^T/e is the thermal energy per 
charge. The factor £(AT) takes into account competi- 
tion between holes and electrons. 

An important parameter in Eq, (5) is the inverse 
Debye screen length K 0 , 

Ki^^/iee^T) (6) 

which depends on the effective density of photore- 
fractive charge N cff and the dc dielectric constant 
e& 0 along the direction of the grating wave vector 
K. 

Combining Eqs. (l)-(6), the gain coefficient F 
can be written in the form: 



( 2 ) r=[ A sin 0/(1 +B- 2 sin 2 c?)](cos2^ I /cos0 I ) 



(7) 

where 0 is the external half angle and 0 l is the 
internal half angle between the two incident laser 
beams. Here we have assumed that the hole-electron 
competition factor f (K) in constant with K in 
order to simplify the data analysis. 

The parameter A is proportional to the effective 
Pockels coefficient r 



erf' 



a = r eff a k ) eXl B = dr/ze I fl= o 



(8) 
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Fig. 3. Experimental curves of T vs. 2$ where the input beams 
wavelength is 488.0 nm. (a) No. 8 sample of CeO 2 (0J0 
niol,%),Fe : Oj(0,05 mol.%):LiNb0 3 ; (b) No. 5 sample of 
CeO 2 (0.10 mol.%XFe 2 O 3 (0.02 mol.%):LiNb0 3 with reduction 
treatment. 



the effective charge density jV eff . The transmitted 
powers of both the weak signal beam and the strong 
pump beam were measured with and without cou- 
pling. The beam diameters were sufficiently large of 
2.5 mm so that the interaction length L of the two 
beams was limited by the physical length of the 
sample. The pump/ signal intensity ratio / p // s was 
3260, the pump beam intensity / p is 4.62 W/crn 2 , 
The input beams was extraordinary polarization. This 
range was a compromise that would minimize pump 
depletion. An iris with a diameter of 1.0 mm was 
placed in front of the power detector to minimize the 
influence of scattering light and to stop all higher 
diffraction lights. 



and is determined by the slope of the plot of r vs. 
29 near 0 = 0. The parameter B is related to the 
effective photorefractive charge density N c{{ : 



B = 



Ait 



4tt V 



N m 



eff 



sin 0 r 



peak 



(9) 



and is determined by 20 peak , the crossing angle at 
which the gain r (20) reaches its maximum value. 
We can obtain the effective photorefractive charge 
density N e{f , an important parameter, from the exper- 
imental curve of r vs. 29: 



MS)' 



w 0 ft B 7 8in^ Mak . 



(10) 



5. Two-beam coupling: experimental results and 
discussion 

The experimental set-up is shown in Fig* 2. An 
Ar + ion laser and a He-Ne laser were used as light 
sources respectively. The two-beam coupling gain 
coefficient T was measured in Ce:LiNb0 3 , 
Fe:LiNt>0 3 , and Ce,Fe:LiNb0 3 crystals to determine 



5.1. Dopant dependence 

Fig. 3 shows the dependence of the steady-state 
two-beam coupling gain coefficient r on the cross- 
ing angle between the two beams 29, They are not 
compared well with the theoretical plot. It is due to 
the stronger light crawling effect in the thinner sam- 
ples [10]. So the crystals show a large angle range in 
which the gain coefficient is higher and the plots 
look flatter. The parameter N eff can be determined 
from Eq. (10) by using some related parameters of 
doped LiNb0 3 crystals as follows: e = 1.602 X 
10" 19 C, 7*= 300 K, A =4.88 (or 5.145, or 6328) 
X 10" 5 cm, k B = 1,38 X 10~ 23 J/K, and ee Q = 
2.832 X 10" 12 C 2 /J cm. N cf{ values determined for 
double and singly doped crystals are listed in Table 
3. 

Compared with the Ce:LiNb0 3 , the N cff value of 
CeJ ? e:LiNb0 3 dramatically increases. While com- 
pared with that of Fe:LiNb0 3 , the increase of that of 
Ce*Fe:LiNb0 3 is not so obvious- So Fe ions still as 
primary photorefractive centers play important roles 
in crystal, and Ce ions as assistants play roles. 
Overall, the photorefractive effect of Ce,Fe:LiNb0 3 



Table 3 

Photorefractive parameters of Ce,Fe:LiNb0 3 crystals 



No. 


Crystal 


Treatment 


Wavelength A(nm) 


Absorption a (cm -1 ) 


Gain JT (cm - 1 ) 


26 (°) 




2 


Ce:LiNb0 3 


reduction 


488.0 


3.6 


24.0 


9.0 


0.2 X 10 15 


3 


Fe;LiNb0 3 


reduction 


488.0 


21.1 


36.4 


31.0 


2.2 X 10 IS 


S 


Ce,Fe:LiNb0 3 


reduction 


488.0 


21.2 


41.7 


36.0 


2-9 X 10 15 
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Table 4 

Photorefractive parameters of Ce,Fe:LiNbO ? crystals 



No. 


Treatment 


Wavelength A (nm) 


Absorption a (cm J ) 


Gain r (cm J ) 


2$ n 


A/ cff (cm 5 ) 


5 


reduction 


488.0 


15.7 


28.4 


28.4 


1.8 X 10 15 


6 


as grown 


488.0 


n.3 


22.2 


20.6 


1.0 x 10 15 


7 


oxidation 


488.0 


1.6 


9.8 


18.8 


0.8 X 10 JS 


8 


reduction 


488.0 


21.2 


41.7 


36.0 


2.9 x 10 15 


9 


heavy reduction 


488.0 


39.0 


38.6 


44.8 


4.4 X 10 15 


8 


reduction 


514.5 


16.2 


39.7 


30.5 


!.9 X 10 15 


8 


reduction 


632.8 


0.7 


27.8 


8.7 


0.1 X 10 15 



was enhanced compared with that of singly doped 



where and <r h are photoionization cross sections 
-of— electron-^and- hole ,- respective! y .- 1 f - there - is only - 
one kind of charge carrier, Eq. (12) can be simplified 

as 

a~<r(Nn+N A ). (13) 

The effect of absorption on N e{{ is due to the concen- 
tration difference between donor and trap (N u 4- N A ). 
As the wavelength becomes longer* the absorption of 
crystal is lower and lower, so, the value of 
decreases (Table 3). 

6. Conclusion 

In summary, we have grown a series of LiNb0 3 
crystals co-doped with Ce- and Fe-ions in different 
levels and fabricated them into experimental samples 
available to photorefractive property measurement. 
The absorption spectra were studied in the wave- 
length range of 300-700 nrn. Compared with the 
singly doped LiNb0 3 , the photorefractive effect of 
Ce,Fe:LiNb0 3 was enhanced. With higher dopant 
concentration in the crystal, its absorption and two- 
beam coupling gain coefficients increase. Oxidation 
and/or reduction treatments have a significant effect 
on photorefractive properties of crystal. To get the 
optimal photorefractive properties, both the dopant 
concentration and proper oxidation/reduction condi- 
tion must be taken into consideration. Double doped 
CcFe:LiNb0 3 crystal shows itself a very promising 
photorefractive material. 
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Table 4 shows the effect of dopant concentration 
on the photorefractive properties of Ce,Fe:LiNb0 3 
crystal. Comparing the experimental results of No, 5 
and No. 8 samples, we can find that the larger the 
dopant concentration, the larger the 7v" cff value. 

5.2. Oxidation / reduction treatment dependence 

II is obviously that A^ ff values are affected by the 
oxidation/reduction treatments greatly. The effec- 
tive charge density N cft can be written as [18] 

K«- • (ID 

Ii actually is the relative concentrations of the donor 
center N D and trap center N A , Total dopant concen- 
tration N — N D +N A . So, proper oxidation/reduc- 
lion treatment will change the photorefractive effi- 
ciency more effectively than the variety of the dopant 
concentrations, because the proper donors and traps 
are the key important factors in photorefractive pro- 
cess. And the photo-charge in crystal transports pri- 
marily by the diffusion when the concentration of Fe 
dopant is lower and the ratio of Fe 2+ to Fe 3+ is 
higher; While the photo-charge transports primarily 
by the photovoltaic effect when the concentration of 
Fe dopant is higher and the ratio of Fe 2+ to Fe 3+ is 
lower [19], 

5.3. Absorption dependence 

Photorefractive absorption coefficient a is de- 
fined as [18] 

* = o c N D -<r h N A (12) 
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